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1.  INTRODUCTION 


The  "Optical  Detection  System  Model"  (ODS  Computer  Code)  was  develop¬ 
ed  under  Contract  F30602- 80-C-0120  to  Rome  Air  Development  Center  (RADC). 
The  overall  effort  was  intended  to  provide  analysts  with  a  computerized 
algorithm  to  simulate  and  evaluate  the  temporal  and  wavelength  response  of 
electro-optical  detection  systems. 

1.1  Overview  of  Optical  Detection  System  Model  Computer  Code 

Ideally,  a  computerized  electro-optical  detection  system  should  include 
every  important  phenomenon  encountered  by  the  radiation  as  it  proceeds  from 
the  source  to  the  sensor,  is  processed  by  the  sensor,  recorded,  and  displayed 
on  a  graphics  terminal.  A  sample  scenario  is  shown  in  Figure  1.1  for  the  de¬ 
tection  of  a  pulse  from  a  laser  as  it  is  scattered  and  attenuated  by  the  atmos¬ 
phere  . 

The  development  of  a  computer  code  for  the  simulation  of  electro-optical 
detection  systems  must,  at  a  minimum,  meet  three  basic  requirements: 

o  The  computer  code  must  be  user  oriented 
o  The  computer  code  must  be  of  modular  construction,  and 
o  The  computer  code  must  be  computationally  efficient. 

The  ODS  computer  code  was  designed  around  these  basic  requirements. 

1.1.1  User  Oriented 

The  key  to  any  successful  computer  code  is  that  it  be  easy  to  use  both 
by  the  individual  generally  unfamiliar  with  the  code  and  those  who  routinely 
use  the  code.  This  requires  that  the  input  data  be  easily  understood  and 
concise.  This  aspect  of  the  ODS  code  development  was  given  the  primary 
priority.  The  user  is  presented  with  a  series  of  "PROMPTS"  on  the  screen  of 
the  CRT  terminal.  Simply  by  answering  a  series  of  questions  the  user  is  led 
through  such  items  as  the  optical  filter  selection,  detector  material  selection, 
amplifier  type,  frequency  response,  etc.  In  addition,  a  number  of  "standard" 
input  signals  can  be  generated  by  the  ODS  code  providing  both  the  experienced 
and  inexperienced  user  with  the  selection  of  an  arbitrary  source  signal  (to  be 


input  by  the  user)  or  the  use  of  one  of  the  supplied  source  signals.  The 
supplied  signals  are  such  that  the  intensity  and  time  may  be  scaled  by  the 
user  to  provide  a  wide  range  of  possible  signal  temparal/intenaity  histories. 

The  inexperienced  user  (or  even  the  experienced  user)  can  ask  for 
assistance  from  the  computer  at  each  entry  point  into  many  of  the  system 
modules.  A  message  appears  on  the  CRT  screen  which  has  information  perti¬ 
nent  to  the  input  being  considered.  The  experienced  user  can  either  employ 
the  "PROMPT"  approach  or  can  employ  input  data  files  already  set  up  through 
previous  use  of  the  ODS  code.  Of  course,  the  experienced  user  may  have 
the  option  of  using  any  input  file  at  his  disposal  for  running  the  code. 


1.1.2  Modular  Construction 

The  computer  code  is  constructed  in  modules  which  represent  different 
facets  of  the  system.  In  this  way  as  more  options  are  developed  these  can  be 
inserted  into  the  code  with  a  minimum  of  difficulty.  An  excellent  example  of 
this  feature  is  the  selection  of  the  detector  material.  A  wide  range  of  detec¬ 
tor  materials  has  been  incorporated  as  part  of  the  resident  data  file  in  the 
ODS  code.  However,  as  new  detector  materials  are  desired  these  new  materials 
can  be  incorporated  into  the  code  with  relative  ease. 

The  modules  currently  employed  include: 

o  SOURCE  MODULE  Permits  the  user  to  select  any  one  of 

a  number  of  stored  source  signals 
such  as  sine  waves,  square  waves, 
step  input ,  etc . ,  or  it  permits  the 
user  to  input  his  own  signal. 

Permits  the  user  to  specify  atmos¬ 
pheric  conditions  for  modifying  the 
signal.  This  routine  currently 
copies  the  input  signal  to  the  output. 

It  is  only  provided  as  a  means  of 
program  extension  at  a  later  date. 

o  OPTICS  MODULE  Permits  the  user  to  specify  the  type 

of  optical  material  used  in  the  collector 
to  determine  the  attenuation  of  the 
signal  at  each  wavelength.  This 
module  also  contains  spectral  filter 
and  etalon  characteristics. 


o  ATMOSPHERIC 
MODULE 
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Permits  the  user  to  choose  one  of  three 
modulating  techniques:  chopper,  circular 
variable  filter,  or  Fourier  Transform 
Spectroscopy. 

Permits  the  user  a  choice  of  visible  and 
infrared  detector  materials  which  includes 
spectral  response  as  well  as  frequency 
response. 

Permits  the  user  a  selection  of  amplifier 
response  functions  including  low  pass, 
high  pass  and  band  pass  amplifiers. 

In  addition  to  the  modules  indicated  above,  each  of  which  can  be  run  as 
a  separate  computer  program,  a  set  of  UTILITY  programs  are  available  which 
enables  the  user  to  EDIT  files  and  PLOT  data.  All  of  the  basic  modules  will 
be  discussed  in  the  following  sections.  The  UTILITY  programs  are  discussed 
in  detail  in  the  Program  Users  Manual  (Ref.  1)  and  the  Program  Maintenance 
Manual  (Ref.  2). 

1 . 2  Application 

It  is  anticipated  that  the  primary  application  of  the  computer  code  is  the 
simulation  of  the  temporal  response  of  sensor  to  an  arbitrary,  time  dependent 
input  signal.  For  this  reason  special  consideration  has  been  given  to  the  fre¬ 
quent  use  of  graphs  of  the  signal  before  it  enters  a  particular  module  and  as 
it  exits.  In  this  way  the  influence  of  the  particular  module  on  the  temporal 
characteristics  of  the  signal  (i.e.,  distribution)  can  be  readily  observed.  In 
addition,  the  construction  of  the  ODS  code  greatly  simplifies  the  selection  of 
preamp /amplifier  design,  for  instance,  since  the  user  may  examine  the  output, 
change  the  design,  and  rerun  only  that  part  of  the  code  in  a  matter  of  min¬ 
utes. 

It  is  further  anticipated  that  the  code  may  be  used  to  simulate  electro- 
optical  data  in  the  presence  of  noise  to  examine  data  processing  techniques. 
Therefore,  care  was  taken  to  insure  that  the  frequency /phase  information  con¬ 
tent  of  the  signal  was  unchanged.  For  this  reason,  the  impulse  invariant 
technique  was  used  for  treating  the  preamp /amplifier  impulse  response  func¬ 
tions.  The  limitations  of  this  technique  are  discussed  in  Section  4.5. 

The  ODS  code  provides  a  unique  method  for  introducing  "new”  engineers 


o  MODULATOR 
MODULE 


o  DETECTOR 
MODULE 


o  AMPLIFIER/ 
PREAMPLIFIER 
MODULES 
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to  the  temporal  response  characteristics  of  electro-optical  sensors.  In  addition, 
with  several  modifications,  the  ODS  code  could  be  employed  to  assist  in  the 
overall  design  of  electro- optical  sensors.  These  modifications  are  discussed 
in  Section  5. 
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2.  SYSTEM  SUMMARY 


This  section  provides  a  brief  summary  of  the  overall  ODS  computer  code 
with  respect  to  systems  requirements,  data  bases,  and  input /output  require¬ 
ments.  Additional  information  may  be  found  in  the  Program  Maintenance 
Manual  and  Program  Users  Manual. 

2. 1  System  Application 

The  Optical  Detection  System  (ODS)  Model  computer  code  is  designed  to 
be  used  interactively  by  the  operator  in  the  evaluation  of  overall  performance 
of  electro -optical  sensors.  The  ODS  code  is  designed  such  that  the  operator 
can  analyze  any  portion  of  the  electro- optical  sensor  such  as  the  filters,  am¬ 
plifier  characteristics,  detector  characteristics,  and  sensor  modulation  charac¬ 
teristics.  The  operator  may  either  analyze  each  component  independently,  as 
groups  of  components,  or  as  the  entire  electro-optical  sensor.  The  operator 
may  employ  stored  signals  for  testing  each  component  or  may  provide  his  own 
particular  signal. 

2. 2  System  Organization 

The  ODS  computer  code  consists  of  20  different  tasks.  Each  of  the 
modules  mentioned  earlier  can  be  run  separately,  as  well  as  programs  for  con¬ 
volution,  editing,  display,  plotting,  and  summing  output  for  multiple  wave¬ 
lengths.  The  entire  control  of  the  overall  operation  of  the  ODS  code  is  main¬ 
tained  through  an  indirect  command  file  -  ODSRUN.  The  execution  of  the 
command  OODSRUN  initiates  a  series  of  commands  which  enables  the  operator 
to  establish  any  combination  of  running  modules  including  the  generation  of 
i  plot  files.  The  first  pass  through  the  command  file  enables  the  user  to  pro¬ 

vide  the  necessary  input  information  required  by  each  of  the  modules.  After 
the  INPUT  phase  the  operator  may  edit  any  of  the  data  files  through  the  EDIT 
command  in  the  ODSRUN  command  file.  When  the  user  is  satisfied  with  the 
input  he  can  execute  the  overall  computation  process  through  the  RUN  command. 

The  organization  of  the  ODS  code  around  a  master  command  file  provides 
great  flexibility  in  adding  new  modules  and  in  placing  the  code  on  other  com¬ 
puters.  The  primary  changes  required  to  place  the  code  on  the  UNIVAC,  for 
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instance,  are  made  in  the  ODSRUN.CMD  command  file.  In  this  file,  the 
appropriate  statements  pertinent  to  the  UNIVAC  can  be  inserted  and  the  same 
overall  structure  of  the  code  employed.  The  different  programs  accessed  by 
the  ODSRUN  command  file  are  indicated  in  Table  1.  The  ODS  code  also  con¬ 
tains  three  data  bases,  one  each  for  the  OPTICS  module,  DETECTOR  module, 
and  AMPLIFIER /PREAMPLIFIER  modules.  These  data  bases  can  be  easily  amend¬ 
ed  to  contain  additional  data. 


TABLE  1 

Running  Tasks  Accessible  to  ODSRUN  Command  File 


PROGRAM 

FUNCTION 

ODSWV 

Runs  SOURCE  module 

WVINP 

Provides  input  to  SOURCE 

ODSATM 

Runs  ATMOSPHERE  module 

ATMINP 

Provides  input  to  ATMOSPHERE 

ODSLNS 

Runs  OPTICS  module 

LNSINP 

Provides  input  to  OPTICS 

ODSMOD 

Runs  MODULATOR  module 

MODINP 

Provides  input  to  MODULATOR 

ODSDET 

Runs  DETECTOR  module 

DETINP 

Provides  input  to  DETECTOR 

ODS H VP 

Runs  PREAMPLIFIER  module 

HVPINP 

Provides  input  to  PREAMPLIFIER 

ODSHVY 

Runs  AMPLIFIER  module 

HVYINP 

Provides  input  to  AMPLIFIER 

ODSEDR 

Runs  EDIT  routine 

ODSDSP 

Runs  DISPLAY  routine 

ODSTOT 

Runs 

routine  which  sums  results  for  multiwavelength 

CNVL 

Runs 

CONVOLUTION  routine 

ODSGRF 

Runs 

PLOT  routine  which  sets  up  plot  file 

ODSXQT 

Runs 

main  module  for  selecting  sequence 

2. 3  Data  Base 

The  ODS  computer  code  employs  three  data  bases.  The  OPTICS  module 
requires  a  data  base  which  contains  the  basic  transmission  characteristics  of 
the  various  optical  materials  and  filters.  A  sample  portion  of  the  data  set  is 
given  in  Table  2.  The  first  line  in  each  data  set  contains  the  code  name  of  the 
material,  the  thickness  (mm),  the  lower  wavelength  (microns),  upper  wavelength 
(microns),  and  the  number  of  wavelength -transmissivity  data  pairs.  The  next 
lines  are  the  wavelength-transmissivity  data  pairs.  In  all  cases  the  wavelength 
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TABLE  2 

SAMPLE  CONTENTS  OF  OPTICAL  MATERIAL  FILE  ODSOPT.INP 
SRF2  3.000  0.120  13.000  9.000 


VYC7905 


INFRASIL 


SAPPHIRE 


3.000 

0.120 

0.120 

0.100 

0.130 

0.600 

0.200 

0.850 

0.300 

0.910 

0.630 

0.940 

0.850 

0.950 

10.000 

0.900 

12.000 

0.500 

13.000 

0.030 

2.000 

1.000 

1.000 

0.925 

2.000 

0.925 

2.600 

0.925 

2.800 

0.930 

3.000 

0.920 

3.200 

0.850 

3.400 

0.550 

3.600 

0.100 

3.800 

0.230 

3.950 

0,100 

4.150 

0.200 

4.300 

0.090 

4.400 

0.100 

4.500 

0.050 

4.700 

0.000 

10.000 

1.000 

1.000 

0.920 

2.000 

0.920 

2.600 

0.920 

2.700 

0.810 

2.800 

0.850 

3.000 

0.840 

3.100 

0.830 

3.300 

0.810 

3.600 

0.630 

3.900 

0.200 

4.000 

0.170 

4.300 

0.000 

1.000 

0.200 

0.200 

0.640 

0.230 

0.770 

0.330 

0.820 

0.700 

0.850 

3.300 

0.850 

4.000 

0.870 

6.000 

0.700 

6.300 

0.100 

4.700  IS. 000 


4.300  12.000 


6.300  8.000 
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is  in  microns. 

The  DETECTOR  module  requires  a  data  base  which  contains  the  properties 
of  the  various  detector  materials  currently  available.  A  sample  of  the  type  of 
data  may  be  seen  in  Table  3.  The  first  line  of  each  data  set  in  this  file  con¬ 
tains  the  detector  temperature  (°K),  time  constant  (seconds),  and  the  number 
of  wavelength-D*  data  pairs.  The  next  lines  are  the  wavelength  (in  microns) 

-  D*  data  pairs.  D*  is  the  detectivity  of  the  detector. 

The  AMPLIFIER  and  PRE -AMPLIFIER  modules  both  access  a  data  file  which 
contains  the  information  necessary  to  generate  a  response  function  for  an  amp/ 
preamp  using  the  Heaviside  routines  as  indicated  in  Table  4.  The  first  line 
of  the  file  contains  the  code  name  of  the  amplifier  /preamp,  and  the  number  of 
poles,  number  of  zeroes,  multiplicative  constant  in  the  quadratic  term,  pre¬ 
factor,  recommended  maximum  time  duration  of  the  response  function  and  the 
recommended  minimum  time  step.  The  time  parameters  are  in  seconds.  The 
next  lines  are  the  poles  and  then  the  zeroes  of  the  amplifier /preamp. 

2.4  General  Description  of  Inputs,  Processing,  Outputs 

The  inputs  required  for  the  ODS  code  are  supplied  either  by  the  user  in 
an  interactive  manner  or  through  the  use  of  data  stored  in  table  form . 

Although  the  information  required  by  each  module  must  be  initially  input  by 
the  user,  the  stored  input  file  may  be  used  for  additional  computer  runs.  In 
the  following  paragraphs  we  will  discuss  in  detail  the  input  required  by  each 
module . 

2.4.1  SOURCE  MODULE 

To  operate  the  ODS  code  the  user  must  specify  the  following  information 
whether  through  prompts  or  by  a  previously  prepared  input  file 


SIGNAM  - Name  of  source  signal  (see  below) 

TSTEP  - TIME  STEP  (time,  in  seconds,  at  which  signal  is 

sampled).  This  is  not  arbitrary  and  must  be 
coordinated  with  the  requirements  set  by  the 
AMPLIFIER  /PREAMPLIFIER  MODULE 

TMAX  . Duration  of  input  signal  (seconds) 

WAVEL  . Wavelength  of  radiation  (microns) 
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TABLE  3 

SAMPLE  CONTENTS  OF  DETECTOR  FILE  ODSDET.INP 

PBS  7.700E401  5.000E-03  U.000 

S.000E-01  4 . OOOE  4 ] 0 
1 . OOQF.400  A.100E410 
1.500E400  8.800E410 
1.750F.400  1.020F.411 
2.000E400  1.200E411 
2.2ej0E400  1.300E411 
2.500E400  1 .600E411 
2.750E.+00  1.800E411 
3.000E400  2.000E411 
3.250E400  2.200E411 
3.500E400  1.900E411 
3.750F400  1 . 600E4 1 1 
4.000E400  9.800E410 
4.250E400  3.600E410 
4.500E400  1.200E410 
4 • 750E400  4.000E409 

PBS  1.930E402  3.000E-03  13.000 

5.000E-01  1.400E411 
7.500E-01  2.000E+11 
1.000E400  2.600E411 
1.500E400  3.700E411 
2.000E400  5.200E411 
2 • 250E400  6.000E411 
2.500E400  6.700E411 
2.750E400  7.000E+11 
3.000E400  5.900E411 
3 . 250E400  3.700E411 
3.500E400  1.000E411 
3.750E400  2.000E410 
4.000E400  2.800E409 

PBS  2.980E402  3.500E-04  14.000 

5.000E-01  9.000E410 
7.500E-01  1.050E411 
1.000F.400  1.200E411 
1.250E400  1.400E411 
1.500E400  1.800E411 
1.750E400  2.000E411 
2.000E400  2.200F.411 
2.170E400  2.100E411 
2.210E400  2.000E411 
2.250E400  1 .800E411 
2.380E400  1 .300E411 
2.500E400  8.000E410 
2.750E400  1.080E410 
3.000E400  7.000E409 


TABLE  4 

SAMPLE  CONTENTS  OF  AMPLIFIER /PRE-AMPLIFIER  FILE 

ODSAMP.INP 


g  6.000E-I00  4.000EI00  7.230E-01  8.960E+10  5.000E-03  1.250E-04 
7.854E+03 
4.712E+01 
3  •  098E  103 
2.513E103 
2.513E+03 
3.378E+03 
O.OOOE-Ol 
O.OOOE-Ol 
O.OOOE-Ol 
O.OOOE-Ol 

H  6.000E+00  4.000E+00  5.400E-01  2.340E+07  2.000F.-01  5.000E-03 
7.8S4K+03 
2.500E-01 
1.092E+02 
6.333E+01 
4.097E+01 
5.460E+01 
O.OOOE-Ol 
O.OOOE-Ol 
O.OOOE-Ol 
O.OOOE-Ol 

F  6.000E+00  4.000E+00  7.230E-01  2.680E+10  1.000E-02  2.500E-04 
7.854E+03 
4.712E+01 
1 • 688E+03 
8.143E+02 
8. 143E+02 
1.848E+03 
O.OOOE-Ol 
O.OOOE-Ol 
O.OOOE-Ol 
O.OOOE-Ol 
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None  of  this  input  data  has  a  default  value.  All  are  required.  Once  the 
user  has  selected  the  type  of  source  signal  to  be  employed  for  the  computations, 
he  must  then  supply  the  necessary  characteristics  of  the  signal.  The  user  will 
be  asked  to  provide  the  additional  information  by  specifying  the  code  word  for 
a  particular  signal.  The  available  source  signals  include  a  sine  wave,  square 
wave,  step-up  function,  step-down  function,  triangle  wave,  ripple  on  a  square 
wave  pulse,  white  noise,  and  1/f  noise.  Samples  of  these  signals  are  given  in 
Figure  2.1.  In  addition  to  the  stored  signals  the  user  also  may  use  his  own 
signal  as  input. 

2.4.2  OPTICS  MODULE 

The  OPTICS  module  contains  the  basic  spectral  filter  data  for  wavelength 
region  isolation.  OPTICS  MODULE  requires  only  the  selection  of  the  optical 
filter  material  and  then  certain  information  necessary  for  that  filter.  All  filter 
transmittance s  are  stored  in  a  file  in  the  computer  for  a  particular  filter  thick¬ 
ness.  If  a  different  thickness  is  input,  the  program  will  perform  the  necessary 
computation  to  increase  or  decrease  the  transmittance  of  the  desired  filter.  An 
etalon  type  filter  may  be  used  also.  The  necessary  information  is  given  in 
Table  5.  Sample  transmission  curves  are  shown  in  Figure  2.2. 

In  addition  to  the  basic  optical /filter  materials  presented  in  Table  5,  the 
user  may  decide  to  employ  an  etalon  for  his  basic  filter.  The  mathematical 
formulation  for  the  etalon  is  provided  in  Section  3.1.2.  The  input  require¬ 
ments  for  using  the  etalon  option  are  given  below. 

ODS  code  word  =  ETALON 

CENTERLINE  WAVELENGTH  -  Wavelength  of  the  centerline  of  the  etalon 

(if  0.0,  assumed  =  to  wavelength  input  in 
beginning  of  run). 

TRANSMISSIVITY  - Transmissivity  at  the  line  center 

INDEX  OF  REFRACTION - Effective  index  of  refraction  of  etalon 

(Sample  Values  2.  -  5.) 

HALFWIDTH  - Halfwidth  of  the  etalon  function  at  the 

halfpower  point 

ANGLE  - Angle  of  incident  radiation  (0°  =  normal 

incidence  and  10°  maximum) 


AC  *  SECONDS  > 


Squar*  htirr  With  Ripplt 


Tint  (seconds) 

TIRE  (SECONDS) 

Sttp  Down 

Whitt  Not  at 

Tint  (SECONDS) 


TINE  (SCCONOS) 


Figure  2.1  Samples  of  the  Input  Signals  which  can  be  generated 
by  the  SOURCE  MODULE. 
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2.4.3  MODULATOR  MODULE 

The  MODULATOR  MODULE  permits  the  user  to  select  one  of  three  methods 
of  "modulating"  the  incoming  signal.  These  methods  are  given  below  together 
with  the  required  input  data  for  each  method. 


CHOPPER 

ODS  code  word  =  CHOPPER 

FREQUENCY  - Chops  per  second  with  a  50%  duty 

cycle  and  a  sinusoidal  output 

CIRCULAR  VARIABLE  FILTER 
ODS  code  word  =  CVF 


SCAN  PERIOD  - Time,  in  seconds,  for  the  CVF  to 

move  from  the  upper  wavelength  to 
the  lower  wavelength 

HALFWIDTH  - Halfwidth  function  of  the  filter  which 


specifies  the  CVF  bandpass  for  a  given 
wavelength.  It  represents  1/2  the 
distance  between  those  wavelength 
points  at  which  the  filtered  signal  will 
be  attenuated  by  a  factor  of  1/2. 

It  is  expressed  as  a  percentage  of 
the  wavelength. 


LOWER  The  lower  wavelength  limit  of  the 

WAVELENGTH  CVF  in  microns 

UPPER  The  upper  wavelength  limit  of  the 

WAVELENGTH  CVF  in  microns 


FOURIER  TRANSFORM  SPECTROSCOPY 
ODS  code  word  =  FTS 
CHARACTERISTIC 

WAVELENGTH  - Microns,  should  be  reasonably  close 

to  other  wavelengths 


2.4.4  DETECTOR  MODULE 

The  DETECTOR  MODULE  performs  two  functions:  1)  it  provides  a  relative 
output  which  is  a  function  of  the  wavelength  of  the  incoming  radiation  (through 
D*),  and  2)  provides  a  frequency  response  function.  All  of  the  required  in¬ 
formation  is  stored  in  the  computer  and  can  be  selected  through  the  appropriate 
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detector  designation.  Detector  characteristics  currently  stored  in  the  computet 
are  given  in  Table  6.  Figure  2.3  provides  sample  D*  curves. 

TABLE  6 

Characteristics  of  Selected  Detectors 


DETECTOR 

TEMPERATURE 

(Kelvin) 

WAVELENGTH 

(microns) 

TIME  CONSTANT 
(seconds) 

PbS 

77. 

0.5-  4.75 

5. 0  x 

IQ' 3 

PbS 

193. 

0.5-  4.0 

3.0  x 

10" 3 

PbS 

298. 

0.5-  3.0 

3.  5  x 

10"4 

PbSe 

77. 

1.0-  6.5 

8.0  x 

10"5 

PbSe 

193. 

1.0-  7.0 

5. 0  x 

10" 5 

PbSe 

298. 

1.0-  4.75 

2.5  x 

10-6 

Si 

78. 

0.575-  1.075 

1.75  x  10"6 

S< 

298. 

0.4-  1.15 

1.0  x 

io8 

InSb 

77. 

1.0-  5.5 

6.0  x 

10- 7 

HgCdTe 

196. 

2.0-  4.4 

1.1  X 

10" 6 

2.4.5  PREAMP /AMPLIFIER  MODULES 

The  PREAMP /AMPLIFIER  MODULES  permit  the  user  to  select  one  of  four 
supplied  filter  pre-amps  (low  pass,  high  pass,  band  pass,  and  band  rejection), 
a  user  supplied  response  function,  and  the  supplied  data  file  where  the  ampli¬ 
fiers  are  described  by  poles  and  zeroes.  This  module  uses  the  time  step 
selected  in  the  SOURCE  module. 

The  amplifier  modules  require  different  inputs  depending  on  which  ampli¬ 
fier  type  is  chosen.  The  various  input  requirements  are  discussed  in  general 
below . 

For  the  filter  amplifiers  indicated  in  Figure  2.4  only  the  attenuation  of 
filter  (in  dB)  and  two  frequencies,  labelled  to  and  to  ,  are  required.  No  time 


UftUtLCNTH  (rtlCRONS) 


Figure  2.3  Sample  Detectivity  (D*)  curves  contained  in  the 
Data  Base  for  the  DETECTOR  MODULE. 


BAND  PASS  AMPLIFIER  FILTER 


LOW  PASS  AMPLIFIER  FILTER 


0.0  OB 


ATTENUATION  (DU) 


0.0  DR 


ATTENUATION  (DB) 


BAND  REJECTION  AMPLIFIER  PILTER 


NIGH  PASS  AMPLIFIER  FILTER 


0.0  OB 


ATTENUATION  (DB) 


0.0  OB 


ATTENUATION  (DU) 


Figure  2.4  Characteristics  of  Filters  Employed  in 
PREAMP /AMPLIFIER  MODULES. 

parameters  are  necessary  as  the  length  of  the  response  function  is  limited  to 
512  points. 

The  filter  amplifier  response  functions  are  generated  directly,  while  those 
described  by  poles  and  zeroes  require  use  of  the  Heaviside  function.  The 
general  form  of  the  Heaviside  function  is, 

PfS-QjXS-Qg) . (S'QN} 

(S-WjXS-Wg) . (s-wM_p(s2-2£wM+wM", 

where 

N  =  number  of  zeroes 
M  =  number  of  poles  and  M  >  N 
Q  =  zeroes 
W  =  poles 

P  =  a  multiplicative  prefactor 

£  =  the  multiplicative  factor  in  the  quadratic  term 


Note:  If  £  =  1,  then  the  quadratic  term  becomes  a  multiple  pole  term  of  order 
2.  If  there  is  no  quadratic  term,  £  should  have  the  value  of  zero. 

At  the  present  time,  the  maximum  number  of  poles  allowed  is  7,  and,  there¬ 
fore,  the  maximum  number  of  zeroes  is  6.  This  is  a  minor  constraint  on  the  oper¬ 
ation  of  the  ODS  code  and  could  easily  be  modified.  No  poles  of  zero  value  are 
allowed . 

It  should  be  noted  here  that  all  amplifier  response  functions  are  responses 
to  the  unit  impulse  and  that  all  time  parameters  are  expressed  in  seconds.  Any 
user  generated  response  function  must  conform  to  these  limitations. 
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3.  MATHEMATICAL  FORMULATION  EMPLOYED  IN 
SELECTED  MODULES 

In  this  section  we  will  present  in  detail  the  mathematical  formulation  employ¬ 
ed  in  certain  modules.  These  modules  include,  optics,  modulator,  detector,  and 
pre-amp /amplifier  module.  The  source  module  requires  no  explanation. 

3.1  Optics  Module 

The  optics  module  performs  two  functions.  It  computes  the  transmissivity 
of  the  optical  material /filters  or  the  transmissivity  of  an  etalon.  Then,  it  multi¬ 
plies  the  incoming  signal  produced  by  the  source  module  by  the  computed  trans¬ 
missivity  . 

3.1.1  Optical  Filter  Material 

The  data  file  ODSOPT.INP  contains  a  number  of  filter  material  as  explained 
in  Section  2.4.2.  The  transmissivity  for  each  of  these  materials  as  a  function  of 
wavelength  is  based  on  a  certain  characteristic  filter  thickness  as  provided  in  the 
table.  The  transmissivity  at  other  thicknesses  is  computed  by  the  relation 


t2  =  tTl] 


V*i 


where  equals  transmissivity  for  second  thickness,  equals  transmissivity  from 
table,  x2  equals  thickness  of  materials  and  Xj  equals  thickness  of  material  obtain¬ 
ed  from  table.  This  relation  is  a  direct  extension  or  direct  consequence  of  Beer's 
law.  In  developing  this  relation  we  have  provided  no  shift  in  the  transmissivity 
due  to  the  angle  of  the  incoming  radiation. 
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3.1.2  Etalon 


The  etalon  transmission  characteristics  are  computed  on  the  basis  of  the 
equation 


T. 


1  + 


2(A-AJj  )' 


AA 


0.5 


where  t  is  the  centerline  transmissivity  at  wavelength,  A  ,  and  A  A-  -  is  the 
full  width  at  half  maximum  transmissivity.  The  prime  values  represent  the 
shifted  center  of  the  etalon  and  change  in  A  Aq  because  of  the  non-normal  in¬ 
cidence  of  the  incoming  radiation.  The  shifted  values  are  given  by 


[1 


1 


AA 


0.5 


A  A 


0.5 


r,., 


02A 


u2aa 


0.5 


where  0  is  the  displacement  of  the  incoming  radiation  from  normal  incident  and 
y  is  the  effective  index  of  refraction  of  the  material  in  the  filter  (Ref.  3). 
Typical  values  of  y  for  the  ultraviolet,  visible,  and  infrared  are  2.0,  2.35, 
and  5.0  respectively  (Ref.  4-6). 

The  use  of  the  formulation  above  does  not  permit  a  decrease  in  the  center - 
line  or  peak  transmissivity  for  non-normal  incident  radiation.  This  decrease  in 
peak  transmissivity  is  a  function  of  many  factors  including  A  ,  AAn  _,  0  and  y. 
However,  in  limiting  the  angle  from  normal  incidence  to  only  a  few  degrees, 
typically  less  than  5°  the  change  in  peak  transmissivity  is  not  severe. 


3.2  Modulator  Module 

The  modulator  module  performs  three  rather  separate  functions.  The  first 
function  is  that  of  a  simple  chopping  wheel  employed  to  change  an  essentially 
DC  infrared  radiation  input  into  an  AC  signal  for  the  detector  and  amplifier 
modules.  The  circular  variable  filter  portion  of  the  modulator  module  provides 
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a  time  dependent  signal  which  is  based  on  the  wavelength  of  the  incoming 
radiation  and  the  characteristics  of  the  circular  variable  filter.  The  Fourier 
transform  spectroscopy  (FTS)  routine  simulates  the  output  of  an  interferometer 
as  a  function  of  time  for  a  number  of  different  wavelengths. 

3.2.1  Chopper 

The  chopper  function  simply  multiplies  the  incoming  signal  from  either  the 
optics  module  or  the  source  module  by  a  sinusoidal  time  varying  function  which 
has  values  ranging  from  0.0  to  1.0.  The  user  simply  specifies  the  frequency 
of  chopping.  The  output  of  the  chopper  routine  is  given  as 

IQ(t)  =  I.(t)  sin  (27ift) 

where  IQ( t )  and  I.(t)  are  the  output  and  input  intensities,  t  is  time,  and  f 
is  the  chopping  frequency  (Hz). 

3.2.2  Circular  Variable  Filter 

The  circular  variable  filter  routine  processes  a  time  varying  input  signal 
with  a  time  varying  filter  function  going  from  a  lower  wavelength,  X ^ ,  to  an 
upper  wavelength,  in  a  period,  T.  Functionally,  this  may  be  represented 
by 


lo.\™  =  Ii,A(t)  CVF(t>  Xl*  V  AX0.5) 

The  CVF  function  represents  a  many  valued  array  which  is  computed  on  the 
basis  of  the  lower  and  upper  wavelengths  and  the  response  halfwidth  of  the 
CVF.  The  response  halfwidth  is  based  on  a  percentage  of  the  wavelength. 
For  instance,  a  3%  resolution  is  given  as 


A  Aq  g  —  0.03  X 

The  spectral  response  of  the  CVF  is  determined  by  employing  a  Gaussian  pro¬ 
file  of  halfwidth  AA„  ...  Thus, 

U.  D 


CVF  =e<0'694  <X-V2,(4X0.5>2) 

where  A  is  the  wavelength  of  the  incoming  signal  and  Aq  is  the  time  varying 
wavelength  position  of  the  CVF  as  it  moves  linearly  in  time  from  X^  to  Ag. 
Thus, 


X 

o 


(A, 


-  A j)  ( 


N)  +  Aj 


for  0  <  <  1  where  T  is  the  time  required  to  move  from  Xj  to  X^  (the  CVF 

repetition  rate).  Of  course,  for  t  >  T  the  ratio  t/T  is  suitably  adjusted  by 
the  number  of  cycles  of  the  CVF  since  t  =  0.  In  the  equation  this  is  given  by 
N. 


It  should  be  noted  that  the  entire  wavelength  range  of  the  CVF  is  divided 
into  1024  points  in  the  computer.  The  limitations  of  1024  points  is  due  to  the 
fact  that  the  computer  code  was  designed  to  operate  on  a  mini  computer.  If 
the  code  should  become  operational  on  a  much  larger  computer  or  one  with 
virtual  memory,  then  the  number  of  points  which  can  be  used  in  the  CVF  can 
be  expanded  to  any  multiple  of  two.  Since  the  current  limitation  is  1024  points 
the  user  should  be  very  careful  in  selecting  the  filter  halfwidth  so  that  suffi¬ 
cient  spectral  resolution  may  be  obtained  over  the  entire  wavelength  region  of 
the  CVF,  For  instance,  if  the  wavelength  range  of  the  CVF  is  from  two  to 
five  microns,  then  the  wavelength  region  will  be  divided  into  wavelength  inter¬ 
vals  of  0.0029  microns.  Such  a  spectral  interval  generally  is  sufficient  for  a 
CVF.  However,  if  a  very  narrow  filter  bandwidth  is  chosen  for  the  spectral 
resolution,  as  defined  by  the  halfwidth,  then  the  wavelength  span  of  the 
CVF  should  be  reduced  accordingly  to  provide  the  proper  spectral  resolution. 


3.2.3  Fourier  Transform  Spectroscopy  (FTS) 

Figure  3.1  illustrates  the  typical  mirror  arrangement  employed  in  an  inter¬ 
ferometer  of  a  type  generally  used  in  FTS.  The  mirror  moves  back  and 
forth  providing  an  interference  pattern  at  the  detector 
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The  interference  pattern  is  given  by  the  equation 


i0<t>  =  i1(t>  +  i2(t)  +  2  /qr^cos  (0(t)> 


where  Ij  and  I2  are  the  intensities  from  mirrors  and  M2  as  perceived  at  the 
detector  and  0(t)  is  the  phase  difference  due  to  the  difference  in  path  length 
traveled  by  the  two  rays.  In  the  case  where  Ij  =  I2  =  Ij/2 

Io(t)  =  2  I.(t)  cos2  ^£1 

It  should  be  noted  that  L  is  the  incoming  intensity.  In  this  formulation  we  have 
ignored  the  reflectivities  of  the  mirrors  and  the  transmission  of  the  splitter  plate 
and  the  compensating  plate. 

The  phase  difference,  0{t),  is  given  by  the  relation 


0(t)  =  2it 


D(t) 

”1 


where  D(t)  is  the  time  varying  difference  in  paths  that  the  two  rays  travel.  In 
normal  FTS  operation ,  the  user  has  control  over  the  rate  of  change  of  the  mirror 
position  as  well  as  the  distance  over  which  the  mirror  may  move.  In  this  computer 
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code,  allowing  the  user  to  select  these  parameters  could  produce  an  FTS  out¬ 
put  which  is  entirely  meaningless,  due  to  the  limited  array  size  available  in  a 
minicomputer.  Thus,  it  was  decided  that  the  mirror  should  move  in  one  direc¬ 
tion  for  1024  time  steps.  During  this  movement  in  one  direction  a  certain 
characteristic  wavelength,  input  by  the  operator,  will  produce  exactly  64 
fringes  on  the  detector.  In  this  way  the  computer  forces  a  mirror  displacement 
and  time  dependency  which  is  consistent  with  an  FTS  output.  Once  these  para¬ 
meters  have  been  determined  by  the  computer  employing  the  characteristic  wave¬ 
length,  the  other  wavelengths  selected  by  the  user  are  phased  employing  the 
predetermined  FTS  parameters.  The  net  result  is  a  typical  FTS  output  spec¬ 
trum. 


3.3  Detector  Module 

The  detector  module  essentially  performs  two  functions.  It  multiplies  the 
incoming  signal  by  the  D*  appropriate  for  that  wavelength  and  it  provides  a 
certain  characteristic  time  dependency.  The  wavelength  dependent  D*  is  in¬ 
dicated  in  the  table  ODSDET.INP.  The  temporal  dependency  of  the  detector 
is  embodied  in  a  characteristic  reseponse  time,  xR.  The  response  of  the  de¬ 
tector  to  a  unit  impulse  is  given  by  the  relation 


This  general  formulation  for  the  detector  module  was  decided  after  examin¬ 
ing  the  temporal  behavior  and  frequency  dependence  of  a  wide  number  of  detec¬ 
tor  materials.  Although  it  does  not  exactly  describe  the  temporal  behavior  of 
each  individual  material,  it  does  provide  a  rather  uniform  formulation  of  the  time 
dependency  of  each  of  the  detector  materials. 

3.4  Preamp /Amplifier  Module 

The  preamp /amplifier  module  performs  three  distinctly  different  types  of 
operation,  each  of  which  produces  a  response  to  a  unit  impulse  which  then  is 
employed  in  the  convolution  routine.  The  first  type  of  operation  employs  a 
generalized  digital  filter  which  can  provide  low  pass,  high  pass,  band  pass,  and 
band  rejection  response  functions.  These  were  discussed  in  Section  2.4.3. 


The  second  type  of  operation  employs  basic  electrical  circuit  analysis  and  is 
based  on  user  definition  of  the  circuit  in  terms  of  poles  and  zeroes.  The  pole- 
zero  response  function  is  analyzed  employing  Heaviside  techniques.  The  third 
type  of  operation  employs  a  user  supplied  "gain"  curve  or  response  function  in 
the  frequency  domain. 

The  response  determined  in  the  time  domain  is  transformed  to  the  frequen¬ 
cy  domain  by  employing  an  FFT.  The  resulting  response  function  in  the 
frequency  domain  is  convolved  with  the  FFT  of  the  incoming  signal*  The  in¬ 
verse  FFT  of  the  results  produces  the  output  of  the  amplifier  in  the  time  do¬ 
main. 

3.4.1  Digital  Filters 

Digital  filters,  as  employed  in  this  section,  are  constructed  in  the  frequen¬ 
cy  domain.  The  inverse  FFT  of  the  filter  provides  the  response  function  in  the 
time  domain.  Next,  the  complex  FFT  is  taken  off  the  time  domain  response 
function  and  used  in  the  convolution  routine.  These  seemingly  unnecessary 
steps  are  required  to  preserve  the  phase  of  the  incoming  signal. 


3. 4. 1.1  Low  Pass  and  High  Pass  Filters 

The  algorithm  employed  in  the  low  pass  filter  routine  is  quite  simple .  For 
frequencies  less  than  a  characteristic  frequency,  u^,  the  signals  are  passed 
unattenuated.  For  frequencies  greater  than  another  characteristic  frequency, 
oig*  the  signals  are  passed  after  they  have  been  attenuated  an  amount  ATTEN 
given  in  decibels.  For  frequencies  between  toj  and  cog.  the  signals  are  attenu¬ 
ated  by  an  amount  based  on  a  linear  interpolation  between  and  ^ .  High 
pass  filters  operate  in  the  reverse  fashion  to  the  low  pass  filters  described 
above . 

3.4. 1.2  Band  Pass  and  Band  Reject  Filters 

Band  pass  filters  operate  in  a  manner  very  similar  to  the  low  pass  and  high 
pass  filters.  The  band  pass  is  defined  by  the  frequencies  and  ^ •  Between 
these  two  frequencies  signals  are  passed  unattenuated.  Outside  these  frequen 
cies  the  signals  are  attenuated  by  an  amount  ATTEN. 
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3.4.2 _ Heaviside  Operation 

The  response  functions,  R(s),  of  many  amplifiers  can  be  described  by  a 
series  of  poles  and  zeroes,  where  the  number  of  poles  exceeds  the  number  of 
zeroes.  These  poles  and  zeroes  are  often  expressed  as  a  polynomial  of  the  form 


R(s) 


Z(s ) 
P(s) 


where  Z(s)  is  the  polynomial  of  zeroes 
and  P(s)  is  the  polynomial  of  poles 


Response  functions  of  this  type  can  be  analyzed  in  a  straightforward  manner  by 
using  the  Heaviside  theorems  as  described  in  Ref.  #7  ,  which  is  closely  followed 
here. 

For  example,  for  an  amplifier  with  2  zeroes  (z^  =  0,  -  -4)  and  4  pole 

terms  (the  last  being  a  quadratic  and  poles  2  and  3  being  equal),  then  R(s) 
would  look  like: 


(s-z.Hs-z,) 

R(s)  =  - 4-^ - r 

(s-PjMs-pg)  (s  -  £p4s+p^) 

5  is  a  constant  term 

This  function  contains  all  of  the  basic  forms  handled  by  the  Heaviside  theorems: 

o 

single  linear  root  (s-pj),  multiple  linear  root  (s-p2)  ,  and  single  quadratic  root 
(s2*2?p4s+p4). 

The  contribution  to  R(t)  1  [  ^  ]  made  by  a  single  linear  term 

(s-a)  is  given  by 


Z(a)  at 
Q(a) 


where  a  =  the  pole  in  the  single  linear  term  being  considered  and  where  Q(a) 
is  the  product  of  all  the  factors  of  P(s)  except  the  (s-a)  term.  Therefore,  for 
this  example  if 


Z1  =  °*  z2  =  4 


Pl  =  - 

2  P,  =■  -1,  P  =  n  (with  1 

then 

R(s)  = 

s(s+4) 

O 

(s+l)(s+l)(s+2)(sz+2s+2) 

s(s+4) 

(s+l)2(s+2)(s2+2s+2) 

here 

Z(s) 

=  s  (s+4) 

and 

Q(s) 

=  (s+1)2  (s2+2s+2) 

evaluating 

Z(s)  and  P  (s)  at  s  =p^  , 

Z(Pj) 

=  -2  (-2  +  4)  =  -4 

Q(Pt) 

=  (-2  +  1 )  2(  ( -  2) 2  +  2(  -  2) 

Therefore,  the  contribution  of  the  single  linear  pole  to  the  overall  temporal 
response  function  is 


=  -2 


e 


-  2t 


The  contribution  to  R(t)  = 


r-l  ,  Z(s) 


(s  -  a)m  is  given  by 


made  by  the  multiple  linear  term 


,at 


4>  (a) (m_  l) 

TffTTTi - 


(m-2) 
<P  (a) 
(m-2)! 


1! 


<t>'(  a) 
1! 


.m-2  m-1 

THT2)T  +  <Ma)  "TnTTTT  1 

where  <t> ( s)  =  quotient  of  Z(s)  and  all  the  factors  of  P(s)  except  (s-a)m.  Note 
that  ^  denotes  the  (m-l)**1  derivative  of  4>(s)  evaluated  at  s  =  a.  There¬ 
fore,  for 
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R(s) 


_ s(s+4) _ 

(s+2)(s+l)2(s2+2s+2) 


0 (s)  =  S^S+— 9 -  and  since  m  =  2, 

(s+2) (sZ+2s+2) 


4>'(s)  = 


2s+4 _ 

(s+2)(s2+2s+2) 


s(s-*-4) _ 

(s+2) 2(s  2+2s+2) 


s(s+4) ( 2s+2) 
(s+2)(s2+2s+2) 2 


and 

<p  (p2)  =  *3 
<t>'  (P2)  =  5 

and 

R„  (t)  =  e_t  [  5- 3t] 
p2 

The  contribution  to  R(t)  1  made  by  the  single  quadratic  terra 


((s  +  a)2 


where 


2 

+  b  )  is  given  by 

e’at 

(<J>.  cos  bt  +  4>r  sm  bt) 

<J>.  is  the  imaginary  part  of  4>  (-a  +  ib) 
0  is  the  real  part  of  0  (-a  +  ib) 


and  0 (s )  is  the  quotient  of  Z(s)  and  all  factors  in  P(s)  except  the  quadratic 
term.  Therefore,  for  our  example. 


_  s(s+4) 

0  ( s )  -  - k- 

(s+2)(s+l) £ 

and 

(s+a)2  +  b2  =  s2  +  2s  +  2 
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and 

a  =  1 
b  =  1 

*(-W)  =  <^»<-1<444>.  ■■  =  -  1  -  31 

(l+i+2)(-lH+l) 


therefore 


!>  =  1,  4>.  =  3 

Kr  yi 


R  (t)  =  e  (-3  cos  t  +  sin  t) 


Finally,  the  overall  temporal  response  function  is 


R (t)  =  -2e  2t  +  (5-3t)  e  *  +  e  t  (-3  cos  t  +  sin  t) 


The  ODS  code  computes  the  temporal  value  of  this  impulse  response 
function  in  time  increments  selected  by  the  user.  Then  the  FFT  of  the  re¬ 
sponse  function  is  performed  and  all  subsequent  operations  occur  in  the  fre¬ 
quency  domain  (i.e.,  the  FFT  of  the  response  function  is  convolved  with  the 
FFT  of  the  incoming  signal) .  At  the  end  of  the  computations  the  inverse  FFT 
is  performed  transferring  the  signal  back  to  the  time  domain. 

3.4.3  Gain  Curve 

This  method  of  filter  design  is  essentially  a  straightforward  application 
of  the  fact  that  a  filter  transfer  function  and  its  impulse  response  are  Fourier 
Transform  pairs.  In  this  program  the  frequency  response  of  the  desired  filter 
is  specified  as  a  gain  curve  -  i.e.,  in  the  frequency  domain  -  and  the  result¬ 
ing  impulse  response  is  derived  from  this  transfer  function. 

The  procedure  is  implemented  as  follows:  first,  the  gain  curve  is  convert¬ 
ed  from  decibels  to  a  linear  scale  of  normalized  magnitudes  between  0  and  1. 
These  magnitudes  are  then  used  as  a  set  of  weighting  coefficients  on  the  set 
of  complex  exponentials  of  unity  magnitude  between  0  and  2tt.  This  results  in 
a  synthetic  "Fourier  Transform"  which  is  then  transformed  back  into  the  time 
domain  by  the  inverse  FFT.  Finally,  a  Blackman  window  is  applied  to  the 
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time  domain  sequence  in  order  to  reduce  the  possible  effects  of  truncation  of 

the  impulse  response,  and  to  improve  the  stability  of  the  filter  at  peaks  in  the 
original  gain  curve. 
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4.  ODS  APPLICATIONS 


The  primary  purpose  of  the  ODS  computer  is  the  simulation  of  the  tem¬ 
poral  response  of  a  wide  range  of  eleetro-optical  sensors.  Therefore,  consider¬ 
ably  more  effort  was  devoted  to  developing  the  proper  temporal  algorithms  than 
wavelength  sensitive  functions.  However,  in  order  to  evaluate  a  sensor  in  an 
elementary  manner  it  is  necessary  that  some  attention  be  given  to  the  wave¬ 
length  of  the  incoming  radiation.  Thus,  the  computer  code  will  consider  the 
wavelength  dependency  of  certain  aspects  of  electro- optical  sensor  design.  The 
manner  in  which  this  wavelength  dependence  is  considered  was  discussed  in 
Section  3.  In  this  section  we  will  present  several  typical  applications  of  the 
ODS  computer  code.  The  actual  manner  in  which  the  operator  interacts  with 
the  ODS  code  is  illustrated  in  the  Appendix  for  each  example  given  below.  The 
Appendix  contains  the  actual  exchange  between  the  operator  and  the  ODS  code. 
Additional  examples  are  provided  in  the  Program  Users  Manual  (Ref.  1). 

4.1  Test  Case  3  -  Program  User's  Manual 

The  primary  purpose  of  this  example  is  to  illustrate  the  results  obtained 
when  a  sinusoidal  input  signal  is  chopped  and  passed  through  a  low  pass  filter 
amplifier  with  the  response  function  given  by 

R(s)  =  C/(s  +  490.2) 

where  C  is  a  constant.  For  this  application  we  employed  the  following  modules 

o  SOURCE  MODULE  provided  the  sinusoidal  signal  at 
50  Hz  and  at  a  wavelength  of  2.5  microns, 
o  OPTICS  MODULE  provided  the  wavelength  dependency 
of  the  sapphire  window. 

o  MODULATOR  MODULE  provided  the  chopper  operation.  The 
chopper  rotated  at  500  Hz, 

o  AMPLIFIER  MODULE  provided  the  impulse  response  function 
indicated  above.  This  is  one  of  the  responses  contained  in 
the  ODS  data  base. 
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Once  the  amplifier  was  selected  it  is  necessary  that  the  time  step  for  the 
computation  be  provided.  The  selection  of  the  time  step  depends  on  two  fac¬ 
tors.  First,  the  time  step  must  be  sufficiently  small  to  provide  an  adequate 
representation  of  the  impulse  response.  Second,  the  duration  of  the  impulse 
response  in  the  time  domain  must  be  of  sufficient  length  that  all  of  the  necess¬ 
ary  features  of  the  response  function  are  included.  The  final  restriction  is 
that  the  current  ODS  code  employs  512  points  in  time  to  represent  the  response 
function.  Therefore,  the  choice  of  time  step  and  duration  of  the  impulse  re¬ 
sponse  must  be  carefully  selected.  The  ODS  code  includes  in  the  data  base 
for  the  amplifier  functions  a  minimum  duration  and  a  maximum  step  size.  For 
the  amplifier  above  the  maximum  step  size  is  0.001  seconds  and  the  minimum 
duration  is  0.020  seconds.  These  were  chosen  through  a  trial  and  error  pro¬ 
cedure  in  which  a  number  of  simple  input  signals  were  examined  employing  a 
range  of  step  sizes  and  response  function  duration. 

Figure  4.1  illustrates  the  output  from  the  chopper  and  indicates  the  basic 
sine  wave  input  and  the  effects  of  chopping  at  500  Hz.  Figure  4.2  illustrates 
an  expanded  portion  of  the  signal  after  it  has  been  processed  by  the  amplifier. 
It  is  obvious  that  the  amplifier  is  too  "slow"  to  provide  a  good  representation 
of  the  chopper  output  signal.  Figure  4.3  indicates  0.5  seconds  of  processed 
signal . 

In  this  situation  the  operator  would  realize  that  the  combination  of  the 
selected  chopper  speed  and  the  low  pass  amplifier  is  not  a  good  match.  There¬ 
fore,  the  operator  could  either  use  a  "faster"  amplifier  or  a  slower  chopper 
speed . 

4.2  Test  Case  4  -  Program  User's  Manual 

This  example  illustrates  the  use  of  the  ODS  code  to  simulate  a  circular 
variable  filter  which  was  discussed  in  Section  3.2.2.  For  this  computation, 
we  employed  a  positive  step  input  from  the  SOURCE  MODULE  and  specified 
three  wavelengths,  4.0,  4.2,  and  4.4  microns,  with  relative  intensities  of  1.0, 
2.0,  and  0.5  respectively.  To  further  illustrate  how  the  ODS  code  functions 
we  started  and  stopped  each  signal  at  different  times.  Thus,  the  signal  at 
4.0  microns  started  at  0.0  seconds  and  ended  at  0.1  seconds.  The  signal  at 
4.2  microns  started  at  0.01  seconds  and  ended  at  0.1  seconds.  The  third 
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TIME  (second 


signal  at  4.4  microns  started  at  0.02  seconds  and  ended  at  0.1  seconds.  It 
should  be  pointed  out  that  the  duration  of  the  signal  does  not  have  to  be  a 
specific  value:  it  can  be  a  very  large  value.  The  ODS  code  will  work  on  the 
signal  in  1024  point  segments. 

The  circular  variable  filter  was  chosen  to  provide  100  scans  per  second. 

In  addition,  the  wavelength  of  the  filter  extended  from  3.5  microns  to  5.0 
microns.  The  resolution  of  the  filter  was  2.0%  of  the  wavelength.  The  output 
from  the  CVF  section  of  the  ODS  code  is  shown  in  Figure  4.4.  It  should  be 
noted  that  the  CVF  scans  from  high  wavelength  to  low.  The  individual  charac¬ 
ter  of  each  wavelength  is  obvious  in  Figure  4.4. 

The  amplifier  consisted  of  2  poles  and  1  zero  with 

„  _  100.04(s  +  617.28) 

~  (s  +  628. 93) (s  +  500.0) 

The  output  of  the  amplifier  is  shown  in  Figure  4.5  which  indicates  that  the 
frequency  response  is  too  low.  The  output  obviously  is  distorted  producing 
a  substantial  error  in  the  relative  intensities  at  the  three  wavelengths,  al¬ 
though  the  position  (wavelength  location)  of  the  peak  intensities  has  not  been 
altered . 

4.3  Test  Case  5  -  Program  Users  Manual 

This  case  illustrates  the  use  of  the  Fourier  Transform  Spectroscopy  part 
of  the  modulator  module.  Wavelengths  of  4.0,  4.2,  and  4.4  microns  were 
selected.  The  amplitudes,  start  times,  and  signal  duration  were  identical  for 
all  three  signals.  These  signals  were  provided  by  the  SOURCE  MODULE.  A 
lower  wavelength  of  3.9  microns  was  used  to  establish  the  necessary  parameters 
for  the  FTS  as  described  in  Section  3.2.3.  We  also  selected  a  detector  material 
of  INSB.  The  output  from  the  detector  is  shown  in  Figure  4.6.  The  charac¬ 
teristic  pattern  of  the  FTS  output  is  clearly  observed. 

This  time  a  simple  high  pass  (HP)  filter  was  selected  for  the  amplifier. 

The  filter  parameters  were  defined  to  have  a  rolloff  frequency  of  100  Hz  and 
with  an  attenuation  of  20  db  at  50  Hz.  All  frequencies  below  50  Hz  were  atten¬ 
uated  20  db.  The  attenuation  between  50  and  100  Hz  is  assumed  to  be  a  linear 
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function.  The  output  from  the  high  pass  filter  is  illustrated  in  Figure  4.7. 

Here  the  reproduction  of  the  detector  signal  is  quite  good.  However,  the  sig¬ 
nal  has  been  displaced  in  time  due  to  the  nature  of  the  digital  filter  design. 
Normally,  this  is  not  a  problem.  However,  if  it  becomes  a  problem  only  a  minor 
change  is  needed  in  the  computer  code . 

4.4  Test  Case  6  -  Program  Users  Manual 

This  example  is  probably  more  typical  of  how  the  ODS  computer  code  may 
be  employed.  Here,  a  simple  square  wave  is  employed  as  the  basic  input  sig¬ 
nal.  However,  we  have  added  some  noise  to  the  signal.  Then  we  have  employ¬ 
ed  an  amplifier  and  a  low  pass  filter,  which  acts  like  an  integrator. 

Figure  4.8  indicates  the  square  wave  and  the  associated  noise.  The 
preamp  which  we  selected  had  one  pole  at  490.0  and  no  zeroes.  The  output 
from  the  preamp  is  illustrated  in  Figure  4.9  which  clearly  shows  the  "cleaning' 
due  to  the  low  pass  nature  of  the  preamp  response  function.  Finally  we  passed 
the  output  from  the  preamp  through  a  low  pass  filter  with  a  rolloff  frequency  at 
100  Hz  and  with  an  attenuation  of  30  db  for  signals  of  120  Hz  and  higher.  The 
results  of  this  process  are  indicated  in  Figure  4.10  where  the  clean  nature  of 
the  signal  is  quite  obvious. 

4.5  TROUBLESOME  PREAMP /AMPLIFIER  TRANSFORMS 

Although  most  transforms  for  conventional  preamps  and  amplifiers  can  be 
handled  by  the  ODS  computer  code,  certain  types  of  transforms  cannot  be 
adequately  treated  in  the  current  version  of  the  code.  In  Section  2.4.5  we 
indicated  that  the  maximum  number  of  zeroes  and  poles  are  6  and  7  respective¬ 
ly.  This  is  a  very  minor  limitation  and  can  be  changed  quite  easily  (Program 
Maintenance  Manual,  Ref.  2). 

A  particular  class  of  transforms  exist  which  would  require  a  major  modi¬ 
fication  to  the  code  and  possible  operation  on  a  computer  with  a  memory  of 
greater  than  32  k  words.  This  limitation  concerns  the  magnitude  of  the  various 
poles  and  zeroes  in  the  transform  and  is  a  direct  result  of  the  limited  number 
of  points  used  in  the  ODS  code  to  represent  the  time  varying  transform.  As 
an  example  consider  the  transform 
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pass 


TIME  (seconds) 

Figure  4.10  Effects  of  Low  Pass  Filter. 


R(s)  (s  +  13.)  (s  +  13000.) 

It  is  virtually  impossible  to  adequately  treat  this  transform  employing  the  al¬ 
gorithms  currently  used  in  the  ODS  code.  The  transform  in  the  time  domain  is 
approximated  by  up  to  512  time  steps.  With  poles  of  greatly  varying  magnitude 
it  is  very  difficult  to  properly  represent  the  transfer  function. 

To  circumvent  this  potential  problem  we  have  examined  two  possible  so¬ 
lutions.  First,  the  user  may  employ  a  very  small  step  size  and  compute  the 
response  for  just  the  larger  pole.  Then,  employing  a  user  generated  program 
the  temporal  output  can  be  adjusted  to  a  larger  step  size.  This  step  would  be 
large  enough  to  accommodate  the  lower  pole.  A  second,  more  satisfying 
approach  to  this  ODS  code  limitation  is  to  consider  the  entire  electrical  circuit 
rather  than  element  by  element.  In  the  problem  illustrated  above  the  circuit 
was  followed  by  one  with  the  transfer  function  given  by 

R(s)  = '  (s  +  4ot:t 

Multiplying  both  transfer  functions  together  to  produce  the  overall  transfer 
function  of  the  system  we  obtain 

s 

R(s)  =  (s  +  13. )(s  +~T30007) ( s  +  4007) 

We  have  indicated  the  gain  curve  for  the  overall  system  and  that  for  the 
two  pole  portion  of  the  system  in  Figure  4.11.  It  is  quite  obvious  that  the 
addition  of  the  third  pole  at  400  rolls  off  the  high  frequency  response  faster 
than  the  13,000  pole.  We  also  have  indicated  the  gain  curve  if  we  eliminate 
the  troublesome  pole  at  13,000.  Clearly,  the  response  characteristics  do  not 
change  significantly.  Therefore,  we  can  adequately  represent  the  overall  sys¬ 
tem  with  two  poles;  one  at  13  and  one  at  400.  It  will  be  recommended  in  the 
next  section  that  the  ODS  code  be  modified  to  provide  for  an  alternate  method 
of  treating  transfer  functions  with  greatly  separated  poles. 


Figure  4.11  Gain  Curves  obtained  through  the  use  of  Poles  and  Zeroes 


5.  RECOMMENDATIONS 


The  current  version  of  the  ODS  computer  code  provides  the  user  with 
many  approaches  to  examining  the  temporal  response  of  electro- optical  sensors. 
At  the  end  of  the  effort  we  re-examined  the  approach  we  had  taken  during  the 
development  of  the  code  and  the  discussions  we  had  with  the  users. 

In  retrospect  we  would  have  taken  somewhat  different  approaches 
in  treating  certain  aspects  of  the  computer  code.  However,  the  actual  ap¬ 
proaches  taken  did  provide  for  an  adequate  simulation  of  the  temporal  response 
of  electro-optical  sensors.  At  this  point  we  believe  that  several  modifications 
and  extensions  to  the  code  would  greatly  enhance  the  ability  of  the  code  to 
treat  a  much  wider  range  of  electro-optical  sensors  and  aid  in  the  design  of 
new  sensors.  In  this  section  we  will  discuss  each  of  the  recommended  modifi¬ 
cations  and  extensions. 

5. 1  Overview  of  Modifications  and  Extensions 

We  have  listed  below  the  recommended  changes  in  the  current  version  of 
the  ODS  code: 

2 

o  Input  signal  and  background  in  watts/cm  /str 
o  Add  scintillation  and  attenuation  to  Atmospheric  Module 
o  Add  elementary  optical  design  package  to  Optics  Module 
o  Add  a  routine  to  perform  a  signal  to  noise  analysis 
o  Expand  the  number  of  points  employed  in  the  FTS  routine 
o  Add  a  wider  range  of  detector  materials  and  types 
o  Add  the  ability  to  treat  a  much  wider  range  of  amplifier 
types  including  standard  RLC  circuits. 

5.2  Signal  and  Background  Input  in  Absolute  Intensity  Units 

Since  the  current  version  of  the  computer  code  was  directly  primarily  at 
the  understanding  of  the  temporal  response  of  the  sensor  very  little  attention 
was  given  to  the  absolute  intensity.  Hence,  although  the  code  can  produce  a 
signal  including  noise  of  several  types,  the  signal  is  not  in  absolute  intensity 
units.  Therefore,  the  user  cannot  determine  if  the  signal  he  is  trying  to 
simulate  might  be  too  weak  for  the  sensor  to  detect.  This  is  particularly 


important  with  the  detector  module  and  the  amplifier  modules  since  the  noise 
which  may  be  produced  by  these  modules  is  not  added  to  the  signal  as  it  is 
processed . 

This  addition  also  requires  updating  the  Optics  Module  to  provide  for 
entrance  optics  preliminary  design.  Thus,  the  user  would  have  the  ability 
to  increase  the  diameter  of  the  entrance  optics  if  the  signal  is  very  weak. 

5.3  Additions  to  Atmospheric  Module 

An  important  addition  to  the  Atmospheric  Module  would  be  the  ability  to 
approximate  the  effects  of  atmospheric  scintillation  since  it  is  one  of  the  fac¬ 
tors  leading  to  laser  signal  attenuation,  in  addition  to  scattering  by  aerosols 
and  absorption.  A  relatively  simple  model  could  be  used  to  permit  the  user 
to  estimate  the  effects  of  various  levels  of  scintillation  on  the  observed  signal. 

5.4  Additions  to  Optics  Module 

Currently  the  Optics  Module  contains  transmissivities  as  a  function  of 
wavelength  for  a  wide  range  of  materials.  A  simple  Beer's  law  relation  is  used 
to  attenuate  the  incoming  signal.  The  Optics  Module  should  contain  an  elemen¬ 
tary  optics  package  which  would  enable  the  use  to  select  a  single  lens  or  one 
of  several  types  of  reflecting  optics.  The  user  also  could  select  coating  for 
the  optics.  In  addition,  the  routine  would  perform  a  computation  to  determine 
the  "blur"  circle  formed  by  the  converging  rays.  In  this  way  the  user  would 
have  some  information  concerning  the  field  of  view  of  the  optics  and  the  basic 
detector  size  (which  is  required  for  a  signal  to  noise  analysis). 

5.5  Signal  to  Noise  Analysis 

The  addition  of  a  signal  to  noise  analysis  to  the  ODS  code  would  greatly 
expand  the  possible  applications  of  the  code.  This  analysis  must  include  the 
radiation  "noise"  from  the  many  surfaces  inside  the  radiometer,  including  the 
lens  and  filter  surface.  In  addition,  the  analysis  must  consider  the  detector 
noise  which  is  a  function  of  detector  temperature,  D*,  area,  and  frequency 
bandwidth  of  the  system.  Finally,  the  analysis  must  include  the  input  noise 
level  of  the  preamplifier. 
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This  use  of  this  signal  to  noise  analysis  could  provide  preliminary  eval¬ 
uation  criteria  for  sensor  design  or  performance.  Of  course,  the  analysis 
must  include  the  optical  collecting  area  as  described  in  Section  5.4.  The 
background  noise,  usually  an  important  aspect  of  the  signal  to  noise  analysis, 
was  discussed  in  Section  5.2. 

5.6  Expanded  Fourier  Transform  Spectroscopy 

At  the  present  time  the  spectral  resolution  capable  with  the  FTS  routine 
of  the  ODS  code  is  limited  by  the  number  of  points  which  can  be  used  to 
represent  the  motion  of  the  mirror.  This  part  of  the  code  must  be  expanded, 
either  using  virtual  memory  or  a  larger  computer. 

5.7  Additional  Detector  Material  and  Types 

The  ODS  code  assumes  that  the  output  from  the  detector  module  is  direct¬ 
ly  proportional  to  the  magnitude  of  the  input  signal  to  that  module.  Such  a 
relationship  includes  most  detector  types.  However,  detector  types  which  do 
not  have  a  direct  relationship  with  the  intensity  are  not  considered.  For  in¬ 
stance,  certain  detectors  respond  like  a  differentiator  to  a  change  in  intensity 
or  polarization.  The  ability  to  treat  this  type  of  detector  should  be  incorporat¬ 
ed  into  the  ODS  code. 

5.8  Wider  Range  of  Preamp /Amplifier  Types 

In  developing  the  ODS  code  we  assumed  that  all  amplifier  responses  would 
be  input  in  the  form  of  responses  of  the  system  to  a  unit  impulse.  It  has 
become  apparent  that  many  users  would  like  the  option  to  employ  a  response 
function  based  on  standard  RLC  circuits.  In  addition,  a  wider  selection  of 
amplifier  response  functions  should  be  incorporated  into  the  ODS  code  data 
base  providing  the  user  greater  flexibility  in  simulating  an  electro-optical 
sensor. 

5.9  Alternate  Techniques  for  Treating  Transfer  Functions 

The  problem  discussed  in  Section  4.5  concerning  troublesome  transfer 
functions  could  be  significantly  relieved  if  alternate  techniques  for  treating 
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transfer  functions  were  available  to  the  user.  The  current  version  of  the 
ODS  code  employs  an  impulse  invariant  technique  for  handling  the  transfer 
functions.  This  technique  was  selected  after  discussion  with  potential 
users  since  it  preservs  frequency  and  phase  information.  The  use  of 
many  digital  signal  processing  techniques  requires  that  the  phase  information 
be  retained. 

The  bilinear  transformation  technique  for  treating  transfer  functions 
would  enable  the  user  to  employ  transfer  functions  with  poles  which  are  much 
further  apart  than  does  the  impulse  invariant  techniques.  The  major  difficulty 
with  the  bilinear  transformation  technique  is  that  it  does  not  retain  all  of  the 
phase  information.  For  those  applications  in  which  phase  information  is  not 
important  the  bilinear  transformation  technique  would  be  quite  useful.  The 
inclusion  of  this  technique  in  the  ODS  code  would  require  a  major  rewrite  of 
a  number  of  routines. 
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I'Ol'SRUN 


;  >  warning:  it  is  necessary  to  clear  all  active  r:>s  fills. 

» 

>*  OK  TO  CLEAR  FILES  AUTOMATICALLY7  [Y/N3JY 

>FIP  ODSXQU.CMD F*/DE 

>F‘IP  ODSMVC.SAV/PU 

>F- 1 F*  t .  SAV/NV=t .  OLD !  */RE 

>FIP  *.PRM/NV=*.NUM5*/RE 

>F'IP  ODSPIC.SAV/NV=ODSPIC.»}*/RE 

PIP  —  NO  SUCH  FILE(S) 

DR0:C350»7]0DSPIC.*»I 
>PIP  ODSPAS.NUM=NLJ 
>» 

>» 

>RUN  ODSXGT 


EEC  OPTICAL  DETECTION  SYSTEM.  V  2.0  313 
SCIENTIFIC  TECHNOLOGY  ASSOCIATES  (1931) 

> INPUT.  EDIT,  RUN?  INPUT 
>SELECT  NEW  SEQUENCE?  C Y/N3 .*  Y 


H+  ODS  MODULES  H+ 


SOURCE 

ATHOS 

OPTICS 

MODULATOR 

DETECTOR 

PREAMP 

AMPLIFIER 

PLOT 


ENTER  MODULE  LIST  (TYPE  VEND'  TO  EXIT) 

::  SOURCE  .OPTICS,  MODULATOR ,  PLOT ,  PREAMP ,  PLOT ,  /END 

>HOW  MANY  SOURCE  WAVELENGTHS  (LE  32767)?  1 


A-2 


...THE  MODULES  YOU  HAVE  REQUESTED  ARE! 

SOURCE 

OPTICS 

MODULATOR  (P) 

PREAMP  (p) 


>OK’  ENTER  Y  OR  NT  Y 
SAVE  THIS  SEQUENCE’ 

warming:  this  will  supersede  any  existing  sequence  , 
c y/nj :  r 

-RUN  WVINP 

THIS  IS  THE  UAVE  FORM  SECTION.  DO  YOU  WANT  HELP?  (Y/NK 

N 

ENTER  WAVELENGTH  (IN  MICRONS)  FOR  WHICH  SIGNAL  IS  TO  DE  DETECTED: 
2.5 

ENTER  NAME  OF  SIGNAL  YOU  WISH  TO  USE: 

SINE 

REFORE  SELECTING  THE  TIME  STEP  TO  PE  USED.  PLEASE 
CONSIDER  THE  FOLLOWING: 

AMPLIFIER  !  ALLOWED  TIME  STEP  <DT ) 


W 

! 

DT 

<«  1.25E-04 

SEC 

H 

I 

DT 

<=  5.00E-03 

SEC 

F 

i 

DT 

2.50E-04 

SEC 

C 

i 

DT 

<=  5.00E-04 

SEC 

D 

! 

DT 

<=  1.00E-03 

SEC 

G 

i 

DT 

<=  1.00E-02 

SEC. 

NOW  PLEASE  ENTER  THE  TIME  STEP  (SEC.)t 
0.0005 

NOW  PLEASE  ENTER  THE  UPPER  LIMIT  ON  TIME  DURATION  OF 
THE  SIGNAL  (SEC.): 

0.5 

THIS  IS  THE  SINE  WAVE  MODULE.  PLEASE  ENTER  VALUES  FOR  THE  FOLLOWING: 

FREQUENCY 

50. 

AMPLITUDE 

1. 

THE  LENGTH  OF  TIME  PEFORE  THE  SINE  PEGINS  (DEADZONE) 

0. 

PASELINE 

0. 


A  -  3 


DO  rOU  WANT  TO  ADD  NOISE  TO  THIS  WAVE  (Y/N)? 

N 

/RUN  LNSINP 

THIS  IS  THE  OPTICAL  MATERIAL  SECTION.  DO  YOU  WANT  HELP  (Y/N)7 

N 

THESE  ARE  THE  AVAILABLE  FILTER  TYPES: 


CORNING 

VYCOR 

SILICA 

MGF2 

SRF2 

VYC7905 

INFRASIL 

SAPPHIRE 

MGO 

T 102 

CDSE 

CDS 

CDTE 

SE 

ZNSE 

CAF2 

DIAMOND 

SILICON 

ZNS 

ETALON 

GAAS 

GE 

LIF 

NACL 

YOU  ARE  NOW 

READY  TO 

INPUT  THE 

NAME  OF 

THE  OPTICAL 

MATERIAL  OF  YOUR  CHOICE  <A8>! 

SAPPHIRE 

ENTER  DESIRED  MATERIAL  THICKNESS  (IN  MM) 

(T^f-e  0.0  to  default  to  stored  thicTr.ess)  (ElO.CKO. 

RUN  MODINP 

*****  THIS  IS  THE  MODULATOR  INPUT  SECTION  ***** 

CHOOSE  ONE  OF  THE  FOLLOWING: 

1)  CHOPPER 

2)  CVF 

3)  FTS 

4)  CHOPPER  AND  CVF 

5)  HELP  UlTH  CHOPPER 

6)  HELP  UITH  CVF 

7)  HELP  UITH  FTS 

ENTER  NUMBER  OF  SELECT  ION*.  1 

***  FOR  THE  CHOPPER: 

PLEASE  ENTER  THE  FREQUENCY  OF  ROT AT I0N(HERT2 ) !  500. 

'.RUN  HVPINP 

THIS  IS  THE  PRE-AMP  SECTION.  USING  THE 
HEAVISIDE  OPERATORS.  DO  YOU  WANT  HELP  (Y/N)? 

N 

DO  YOU  WISH  TO  USE  ONE  OF  THE  SUPPLIED  ROUTINES 
(Y/N)? 

Y 

WHICH  RESPONSE  FUNCTION  DO  YOU  WANT?  THESE  ARE  THE  AVAILABLE  CHOICES: 
U 
H 
F 
C 
D 
G 

LP 

BP 

HP 

BR 


A-4 


NOW  ENTER  YOUR  CHOICE: 

Ii 

PLEASE  ENTER  THE  UPPER  LIHIT  ON  THE  TIME  FOR  THE 
IMPULSE’ RESPONSE  (ElO.O).  THIS  UILL  RE  IN  SECONDS. 
SELECT  YOUR  TIME  LIMIT  TO  RE  RETUEEN  THE  FOLLOWING  TWO 
TIMES:  2.560E-01  SECONDS  AND  2.000E-02  SECONDS. 

0.04 

>» 

:*  CONTINUE  WITH  ODS?  CY/N3:Y 
>i 

PUN  ODSXQT 


CCC  OPTICAL  DETECTION  SYSTEM.  V  2.0  33] 
SCIENTIFIC  TECHNOLOGY  ASSOCIATES  (1981) 


>  INPUT.  EDIT.  RUN?  RUN 
SELECT  NEW  SEGUE  NCI''  LY/Nit  N 

..pip  odssum.old/nv=nl: 

>pip  odssum.neu/nv=nl: 

>pip  odsinp. old/nv=nl: 

>PIP  ODS  INF’ .  NEU/NV=NL  • 

>RUN  ODSUU 

>PIP  ODSINP.OLD/NV=OOSINP.NEU.O/RE 
:  P IP  ODS INP . NEU/NU=ODSGUT . NEW i O/RE 
>PIP  ODSWV . NUM/NV=ODSWV . NUM.-l/RE 
>RUN  ODSLNS 

>P IP  ODSINP . OLD/NV=ODSIHP . NEW O/RE 
>PIP  ODSINP .N£U/NV=ODSOUT .NEW, O/RE 
>RUN  ODSMOD 

>PIP  ODSINP . OLD/NV=ODSINP . NEW » O/RE 
>PIP  ODSINP. NEW/NV=ODSOUT. NEW? O/RE 
>RUN  ODSGRr 
>F:UN  ODSTOT 

>PIP  ODSSUM . OLD/NV=ODSSUM .NEW! O/RE 
>PIP  ODSINP , NEW/NV=ODSSUM . OLD 5 O/RE 
>RUN  ODSHVP 

>PIP  ODSIR.NEW=ODSOUT. NEW? O/RE 
>RUN  CNVL 

>P IP  ODSINP . OLD/NV=ODSI NP . NEW ; O/Rf 
>PIP  ODSIR.OLD/NV=ODSIR. NEW; O/RE 
>PIP  ODSINP. NEW/NV=ODSOUT. NEW, O/RE 
RUN  ODSGRF 
i RUN  ODSDSP 
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TEST  CASE  4 
Program  Users  Manual 


A-  7 


POI'SPUN 


.  ;  »  warning:  it  is  necessary  to  clear  all  active  cos  files. 

> » 

>*  OK  TO  CLEAR  FILES  AUTOMATICALLY7  CY/NTiY 

/PIP  ODSXQU.CMD »*/BE 

>F'IF'  ODSMVC .  SAV/PU 

>P1F‘  #.SAV/NV=*.OLDi*/RE 

> F‘ I F*  *.F'RH/NV=*.NUM»*/RE 

>P IF'  ODSPIC . SAV/NV=ODSP I C .  * i */RE 

>pip  odspas.num=nl: 

y\ 

>RUN  ODSXOT 


EEC  OPTICAL  DETECTION  SYSTEM,  V  2.0  ]]] 
SCIENTIFIC  TECHNOLOGY  ASSOCIATES  (1981) 

> INPUT,  EDIT,  RUN?  INPUT 
/SELECT  NEW  SEQUENCE?  E Y/N3 !  Y 

+++  ODS  NODULES  +++ 


SOURCE 

ATMOS 

OPTICS 

MODULATOR 

DETECTOR 

PREAMP 

AMPLIFIER 

PLOT 


ENTER  MODULE  LIST  (TYPE  VEND'  TO  EXIT) 
/SOURCE , hODUL ATOR , PLOT , AMPL IFI ER , PLOT , /END 

:HOU  MANY  SOURCE  WAVELENGTHS  (LE  32767)7  3 


...THE  MODULES  YOU  HAVE  REPUESTED  ARE? 


SOURCE 

MODULATOR  (P) 

AMPLIFIER  (P) 

.OK7  ENTER  Y  OR  N:  Y 

A-8 


i 


:  SAVE  THIS  SEQUENCE7 

:>  warning:  this  will  supersede  any  existing  sequence  . . 

>[ Y/N3 :  Y 
>RUN  WVINP 

THIS  IS  THE  WAVE  FORM  SECTION.  DO  YOU  WANT  HELP?  < Y/N) I 

N 

ENTER  WAVELENGTH  (IN  MICRONS)  FOR  WHICH  SIGNAL  IS  TO  BE  DETECTED 
4. 

ENTER  NAME  OF  SIGNAL  YOU  WISH  TO  USEt 
STEP+ 

BEFORE  SELECTING  THE  TIME  STEP  TO  BE  USED.  PLEASE 
CONSIDER  THE  FOLLOWING: 

AMPLIFIER  !  ALLOWED  TIME  STEP  (DT) 


u 

} 

DT 

<=  1.25E-04 

SEC 

H 

| 

DT 

<=  5.00E-03 

SEC 

F 

| 

DT 

<=  2.50E-04 

SEC 

C 

| 

DT 

<=  5.00E-04 

SEC 

D 

1 

DT 

<=  1.00E-03 

SEC 

G 

| 

DT 

<=  1.0CE-02 

SEC 

NOW  PLEASE  ENTER  THE  TIME  STEP  (SEC.)J 
l.E-4 

NOW  PLEASE  ENTER  THE  UPPER  LIMIT  ON  TIME  DURATION  OF 
THE  SIGNAL  (SEC.)  I 
0.1 

THIS  IS  THE  STEP-PLUS  ROUTINE  (A  STEP  UP). 

PLEASE  ENTER  THE  FOLLOWING  VALUES: 

THE  AMPLITUDE 

1.0 

LENGTH  OF  TIME  BEFORE  THE  SIGNAL  BEGINS  (DEADZONE) 

0. 

THE  BASELINE  (MINIMUM  AMPLITUDE) 

0. 

DO  YOU  WANT  TO  ADD  NOISE  TO  THIS  WAVE(Y/N)? 

N 

>RUN  WVINP 

THIS  IS  THE  WAVE  FORM  SECTION.  DO  YOU  WANT  HELP?  (Y/N): 

N 

ENTER  WAVELENGTH  (IN  MICRONS)  FOR  WHICH  SIGNAL  IS  TO  BE  DETECTED*. 
4.2 

ENTER  NAME  OF  SIGNAL  YOU  WISH  TO  USE*. 

STEP+ 

NOW  PLEASE  ENTER  THE  UPPER  LIMIT  ON  TIME  DURATION  OF 
THE  SIGNAL  (SEC.): 

0.09 
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THIS  IS  THE  STEP-PLUS  ROUTINE  (A  STEP  UP). 

PLEASE  ENTER  THE  FOLLOUING  VALUES: 

THE  AMPLITUDE 

2. 

LENGTH  OF  TIME  BEF  ORE  THE  SIG'iAL  IEGINS  ( DEADZONL  > 

O.Ol 

THE  BASELINE  (MINIMUM  AMPLITUDE) 

0. 

DO  YOU  WANT  TO  ADD  NOISE  TO  THIS  UAVE(Y/N>? 

N 

>RUN  WVINP 

THIS  IS  THE  WAVE  FORM  SECTION.  DO  YOU  WANT  HELP?  ( Y/N) ! 

N 

ENTER  WAVELENGTH  (IN  MICRONS)  FOR  WHICH  SIGNAL  IS  TO  BE  DETECTED 

4.4 

ENTER  NAME  OF  SIGNAL  YOU  WISH  TO  USE: 

STEP+ 

NOW  PLEASE  ENTER  THE  UPPER  LIMIT  ON  TIME  DURATION  OF 
THE  SIGNAL  (SEC.): 

0.08 

THIS  IS  THE  STEP-PLUS  ROUTINE  (A  STEP  UP). 

PLEASE  ENTER  THE  FOLLOWING  VALUES! 

THE  AMPLITUDE 
0.5 

LENGTH  OF  TIME  BEFORE  THE  SIGNAL  BEGINS  (DEADZONE) 


0.02 

THE  BASELINE  (MINIMUM  AMPLITUDE) 

0. 

DO  YOU  WANT  TO  ADD  NOISE  TO  THIS  WAVE(Y/N)? 
N 

>RUN  MODINP 


*****  THIS  IS  THE  MODULATOR  INPUT  SECTION  ***** 

CHOOSE  ONE  OF  THE  FOLLOWING! 

1)  CHOPPER 

2)  CVF 

3)  FTS 

4)  CHOPPER  AND  CVF 

5)  HELP  WITH  CHOPPER 

6)  HELP  WITH  CVF 

7)  HELP  WITH  FTS 

ENTER  NUMBER  OF  SELECTIONS 

U*  FOR  THE  CVF! 

PLEASE  ENTER  THE  SCAN  PERIOD  (IN  SECONDS):  0.01 
ENTER  THE  SPECTRAL  RESOLUTION  IN  7.  OF  UAVEIENGTH1  2.0 
PLEASE  ENTER  THE  LOU  WAVELENGTH!  3.5 
PLEASE  ENTER  THE  HIGH  WAVELENGTH!  5.0 
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, run  hvyinp 

THIS  IS  THE  AMPLIFIER  SECTION.  USING  THE 
HEAVISIDE  OPERATORS.  DO  YOU  UANT  HELP  (Y/N)? 

N 

DO  YOU  WISH  TO  USE  ONE  OF  THE  SUPPLIED  ROUTINES 
(Y/N)? 

Y 

LHICH  RESPONSE  FUNCTION  JO  YOU  UANT?  THESE  ARE  THE  AVAILABLE  CHOICE 

U 

H 

F 

C 

D 

G 

LP 

BP 

HP 

BR 

NOU  ENTER  YOUR  CHOICE: 

C 

PLEASE  ENTER  THE  UPPER  LIMIT  ON  THE  TIME  (IN  SECONDS) 

FOR  THE  IMPULSE  RESPONSE  (ElO.O) 

SELECT  YOUR  TIME  LIMIT  TO  BE  BETWEEN  THE  FOLLOWING  TWO 
TIMES.*  5.120E-02  SECONDS  AND  2.000E-02  SECONDS. 

4.E-2 

>. 

>*  CONTINUE  WITH  ODS?  CY/N3:Y 

>*. 

>RUN  ODSXQT 


ETC  OPTICAL  DETECTION  SYSTEM.  V  2.0  330 
SCIENTIFIC  TECHNOLOGY  ASSOCIATES  (1981) 


> INPUT »  EDIT.  RUN?  RUN 
>SELECT  NEW  SEQUENCE?  CY/N3*.  N 

>pip  odssum.old/nv=nl: 

>pip  odssum.new/nv=nl: 

>F'IP  odsinp.old/nv=nl: 

>pip  odsinp.neu/nv=nl: 

/RUN  ODSWV 

>pip  odsinp. olb/nv=odsinp.neu;o/re 

>P IP  ODSINP . NEW/NV=ODSOUT . NEW i O/RE 
>PIP  ODSWV .-NUM/NV=ODSWV .NUM.-l/RE 
-RUN  ODSMOD 

; P I P  ODS I NP . OL D/NV=ODS I NP . NEW i O/RE 
/PIP  ODSINP . NEW/NV=ODSOUT . NEW » O/RE 
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>RUN  ODSGRF 

:run  ohs tot 

>P IP  ODSSUM . GLD/NV-ODSSUM . NEW f O/RE 
>RUN  ODSWV 

>P IP  ODSINP . OLD/NV=ODS INP . NEW » O/RE 
> PIP  OHS INP . NEW/NV=ODSOUT .NEW? O/RE 
>PIP  ODSUV . NUM/NV=ODSUV . NUM» - 1/RE 

:run  odshoh 

>PIP  OHS INP . OLD/NU=ODSINP .  NEU » O/RE 
>PIP  ODSINP . NEU/NV=ODSOUT .NEU  JO/RE 
>RUN  ODSGRF 
✓  RUN  OHSTOT 

>P IP  ODSSUM . OLD/NV=ODSSUH . NEW i O/RE 
>RUN  ODSUV 

>PIP  ODSINP . OLD/NV=ODSINP . NEW » O/RE 
>PIP  ODSINP . NEU/NV=ODSQUT . NEU » O/RE 
>PIP  ODSWV. NUM/NV=ODSUV.NUM»-l/RE 
>RUN  ODSMOU 

>P IP  ODSINP . OLD/NV=ODS INP . NEU i O/RE 
>PIP  ODSINP. NEU/NV=ODSOUT. NEW! O/RE 
>RUN  ODSGRF 
>RUN  ODSTOT 

>P IP  ODSSUM . OLD/NV=ODSSUM .  NEU ! O/RE 
>PIP  ODSINP. NEU/NV=ODSSUM.OLD50/RE 
>RUN  ODSHVY 

>PIP  ODSIR.NEU=ODSOUT .NEU50/RE 
>RUN  CNVL 

>PIP  ODSINP. OLD/NV=ODSINP, NEU JO/RE 
>PIP  ODS IR.OLD/NV=ODSIR. NEW? O/RE 
>PIP  ODSINP. NEW/NV=ODSOUT . NEU ? O/RE 
>RUN  ODSGRF 
>RUN  ODSDSP 
>5 

>%  CONTINUE  WITH  ODS?  CY/ND1N 
>» 

>?  » 

>»  »  END  OF  ODS  RUN 

>f  » 

>?  »  (YOU  MAY  UISH  TO  PURGE  FILES  *.SAV  AND  *.PRM  AT  THIS  POINT) 

>»  » 

>r 

>(?  <EOF> 

> 

>PIP  *.sav?*/de 

> PIP  *.PF:M?*/DE 

> 
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POI'SRUN 

/» 

>i  :■>  warning:  it  is  necessary  to  clear  all  active  ods  filt:.. 

; » 

>*  OK  TO  CLEAR  FILES  AUTOMATICALLY?  CY/NKY 
>PIP  ODSXQU.CMD ;*/DE 
>F'IP  ODSMVC.SAV/PU 
PIP  —  NO  SUCH  FILE<S) 

SY0:C35O,73ODSMVC.SAV 

>P1P  *.SAV/NV=*.OLD»*/RE 

>PIP  *.PRH/NV=*.NUMi*/RE 

>P1P  ODSPIC.SAV/NV=ODSPIC.*j*/RE 

PIP  —  NO  SUCH  FILE(S) 

DR0:C350»710DSPIC.*»* 

>pip  odspas.num=nl: 

>i 

>t 

>RUN  ODSXQT 


CCC  OPTICAL  DETECTION  SYSTEM,  V  2.0  333 
SCIENTIFIC  TECHNOLOGY  ASSOCIATES  (1981) 


> INPUT »  EDIT,  RUN?  INPUT 
>SELECT  NEW  SEQUENCE?  EY/N3 :  Y 


B+  ODS  MODULES  +++ 


SOURCE 

ATMOS 

OPTICS 

MODULATOR 

DETECTOR 

PREAMP 

AMPLIFIER 

PLOT 


-ENTER  MODULE  LIST  (TYPE  VEND’  TO  EXIT) 

:  SOURCE, MODULATOR, PLOT, DETECTOR, PLOT, AMPLIFIER, PLOT, /END 

>HOW  MANY  SOURCE  WAVELENGTHS  (LE  32767)?  3 


A- 14 


...THE  MODULES  YOU  HAVE  REQUESTED  ARE 


SOURCE 

MODULATOR  (P) 
DETECTOR  (p> 
AMPLIFIER  (P) 


>OK?  ENTER  Y  OR  NJ  Y 
SAVE  THIS  SEQUENCE? 

/i  warning:  this  uill  supersede  any  existing  sequence  : 

>CY/ND :  Y 
>RUN  UVINP 

this  is  the  wave  form  section,  do  YOU  WANT  HELP?  (Y/N>: 

N 

ENTER  WAVELENGTH  (IN  MICRONS)  FOR  WHICH  SIGNAL  IS  TO  BE  DETECTED 
4.0 

ENTER  NAME  OF  SIGNAL  YOU  WISH  TO  USE l 
STEP+ 

BEFORE  SELECTING  THE  TIME  STEP  TO  BE  USED*  PLEASE 
CONSIDER  THE  FOLLOWING: 

AMPLIFIER  !  ALLOWED  TIME  STEP  (DT) 


W 

H 

F 

C 

D 

G 


DT  <=  1.25E-04  SEC. 
DT  <=  5.00E-03  SEC. 
DT  <=  2.50E-04  SEC. 
DT  <=  5.00E-04  SEC. 
DT  <=  1.00E-03  SEC. 
DT  <=  1.00E-02  SEC. 


NOW  PLEASE  ENTER  THE  TIME  STEP  (SEC.) I 

0.001 

NOW  PLEASE' ENTER  THE  UPPER  LIMIT  ON  TIME  DURATION  OF 
THE  SIGNAL  (SEC.)’. 

1.0 

THIS  IS  THE  STEP-PLUS  ROUTINE  (A  STEP  UP). 

PLEASE  ENTER  THE  FOLLOWING  VALUES' 

THE  AMPLITUDE 

1.0 

LENGTH  OF  TIME  BEFORE  THE  SIGNAL  BEGINS  (DEADZONE) 

0. 

THE  BASELINE  (MINIMUM  AMPLITUDE) 

0. 

DO  YOU  WANT  TO  ADD  NOISE  TO  THIS  WAVE(Y/N)? 


/RUN  WVINP 

THIS  IS  THE  DAVE  FORM  SECTION.  DO  YOU  WANT  HELP?  (Y/N>! 

N 

ENTER  WAVELENGTH  < IN  MICRONS)  FOR  WHICH  SIGNAL  IS  TO  PE  DETECTED 
4.2 

ENTER  NAME  OF  SIGNAL  YOU  WISH  TO  USE! 

STEP* 

NOW  PLEASE  ENTER  THE  UPPER  LIMIT  ON  TIME  DURATION  OF 
THE  SIGNAL  (SEC.); 

1.0 

THIS  IS  THE  STEP-PLUS  ROUTINE  (A  STEP  UP). 

PLEASE  ENTER  THE  FOLLOWING  VALUES! 

THE  AMPLITUDE 

2.0 

LENGTH  OF  TIME  PEFORE  THE  SIGNAL  DEGINS  (DEADZONE) 

0. 

THE  BASELINE  (MINIMUM  AMPLITUDE) 

0. 

DO  YOU  WANT  TO  ADD  NOISE  TO  THIS  WAVEU/N)? 

N 


RUN  WV1N  ' 

THIS  IS  THE  WAVE  FORM  SECTION.  DO  YOU  WANT  HELP'  <  Y/W) ! 

N 

ENTER  WAVELENGTH  (IN  MICRONS)  FOR  WHICH  SIGNAL  IS  TO  PE  DETECTED! 

4.4 

ENTER  NAME  OF  SIGNAL  YOU  WISH  TO  USE! 

STEP+ 

NOW  PLEASE  ENTER  THE  UPPER  LIMIT  ON  TIME  DURATION  OF 
THE  SIGNAL  (SEC.)! 

1.0 

THIS  IS  THE  STEP-PLUS  ROUTINE  (A  STEP  UP). 

PLEASE  ENTER  THE  FOLLOWING  VALUES! 

THE  AMPLITUDE 
0.5 

LENGTH  OF  TIME  PEFORE  THE  SIGNAL  PEGINS  (DEADZONE) 

0. 

THE  PASELINE  (MINIMUM  AMPLITUDE) 

0. 

DO  YOU  WANT  TO  ADD  NOISE  TO  THIS  WAVE(Y/N>? 

N 

>RUN  MODINP 

*****  THIS  IS  THE  MODULATOR  INPUT  SECTION  ***** 

CHOOSE  ONE  OF  THE  FOLLOWING! 

1)  CHOPPER 

2)  CVF 

3)  FTS 

4)  CHOPPER  AND  CVF 

5)  HELP  WITH  CHOPPER 

6)  HELP  WITH  CVF 

7)  HELP  WITH  FTS 

ENTER  NUMBER  OF  SELECTIONS 
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n*  FOR  THE  FTSt 

PLEASE  ENTER  THE  LOW  UAVELENGTH!  3.9 
>RUN  DETINP 

THIS  IS  THE  DETECTOR  MODULE.  YOU  WILL  BE  ASKED  TO 
PICK  A  DETECTOR  TYPE  AND  OPERATING  TEMPERATURE 
COMPATABLE  WITH  THE  WAVELENGTH  YOU  HAVE  ALREADY 
CHOSEN.  FAILURE  TO  DO  SO  MAY  BE  CORRECTED  DURING 
THIS  INPUT  SESSION  OR  BY  INVOKING  THE  EDIT  MODI 
AFTER  YOU  HAVE  EXITTED  FROM  THIS  ROUTINE. 

DO  YOU  WANT  HELP?  < Y/N) 

N 

THESE  ARE  YOUR  DETECTOR  CHOICES*. 

PBS 

PBSE 

SI 

INSB 

HGCDTE 

PLEASE  ENTER- ONE  OF  THE  ABOVE  DETECTOR  MNEMONICS! 

INSB 

THESE  ARE  THE  AVAILABLE  TEMPERATURES  FOR  YOUR  DETECTOR 
CHOICE! 

77.000 

PLEASE  ENTER  A  TEMPERATURE  TIN  DEGREES  KELVIN) *. 

77. 

RUN  HVYI f 

THIS  IS  THE  AMPLIFIER  SECTION,  USING  THE 
HEAVISIDE  OPERATORS.  DO  YOU  WANT  I'ELP  (Y/N)7 

N 

DO  YOU  WISH  TO  USE  ONE  OF  THE  SUPPLIED  ROUTINES 
(Y/N)? 

Y 

WHICH  RESPONSE  FUNCTION  DO  YOU  WANT?  THESE  ARE  THE  AVAILABLE  CHOICES! 

U 

H 

F 

C 

D 

G 

LP 

BP 

HP 

BR 

NOW  ENTER  YOUR  CHOICE! 

HP 

IT  WILL  NOT  BE  NECESSARY  TO  ENTER  A  MAXIMUM  TIME. 

IT  WILL  BE  COMPUTED  BASED  ON  A  FIXED  NUMBER  OF  POINTS  =  256. 


NOU  YOU  MAY  ENTER  THE  FIRST  FREQUENCY. 

PLEASE  REMEMBER  THAT  IF  YOU  HAVE  CHOSEN  LP  OR  HP  AMPLIFIERS! 
THAT  THIS  IS  THE  CUTOFF  FREQUENCY.  IF  YOU  HAVE  CHOSEN  THE 
BP  OR  BR  AMPLIFIERS.  THEN  THIS  IS  THE  FREQUENCY  OF  THE 
LEADING  EDGE  OF  THE  BAND.  PLEASE  ENTER  THE  1ST  FREQUENCY: 

100. 

NOU  YOU  MAY  ENTER  THE  2ND  FREQUENCY. 

PLEASE  REMEMBER  THAT  IF  YOU  HAVE  CHOSEN  LP  OR  HP  AMPLIFIERS. 
THAT  THIS  FREQUENCY  UITH  THE  1ST  FREQUENCY  MARK  THE 
TRANSITION  BAND.  FOR  THE  LP.  THIS  UILL  BE  >  1ST  FREQ. 

FOR  THE  HP.  THIS  UILL  BE  <  1ST  FREQ.  IF  YOU  HAVE  CHOSEN 
THE  BP  OR  BR  AMPLIFIERS.  THEN  THIS  IS  THE  FREQUENCY  OF 
THE  TRAILING  EDGE  OF  THE  BAND.  PLEASE  ENTER  THE  2ND  FREQUENCY: 
50. 

PLEASE  ENTER  THE  ATTENUATION  (IN  DB>! 

20. 

>» 

'✓*  CONTINUE  UITH  ODS?  CY/NV.Y 

>; 

>RUN  ODSXQT 


CTC  OPTICAL  DETECTION  SYSTEM.  V  2.0  D31 
SCIENTIFIC  TECHNOLOGY  ASSOCIATES  (1981) 


> INPUT.  EDIT,  RUN?  RUN 
✓SELECT  NEU  SEQUENCE?  CY/N 31  N 

✓pip  odssum.old/nv^nl: 

✓pip  odssum.neu/nv=nl: 

>pip  odsinp. old/nv=nl: 

✓pip  odsinp. neu/nv=nl: 

-RUN  0DSUV 

✓PIP  0DSINP.0LD/NV=0DSINP.HEU;0/RE 
✓PIP  0DSINP.NEU/NV=0DS0UT .NEU50/RE 
✓ PIP  ODSUV . NUM/NV-ODSUV . NUMJ - 1 /RE 
✓RUN  ODSMOD 

>P I P  ODS I NP . 0LD/NV=0DS I NP . NEU i O/RL 
✓PIP  ODSINP .NEU/NV=ODSOUT .NEW? O/RE 
✓RUN  ODSGRF 
✓RUN  ODSDET 

>FIP  GDSIR.NEU=0DS0UT.NEU50/RE 
>RUN  CNVL 

>P IP  ODSINP . 0LD/NV=0DSINP . NEW.O/RE 
>PIP  ODSIR «0LD/NV=0DSIR . NEW. O/RL 
>PIP  ODSINP . NEU/NV=ODSOUT . NEW ? O/RE 
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>RUN  ODSGRF 
>RUN  ODSTOT 

>PIP  ODSSUM . OLD/NV=ODSSUM . NEW » O/RE 
>RUN  ODSWV 

>PIP  ODSINP.OLD/NV=ODSINP. NEW  JO/RE 
/PIP  ODSINP . NEW/NV=ODSOUT . NEW i O/RE 
>PIP  ODSWV . NUM/NV=ODSWV . NUMJ -1/RE 
>RUN  ODSMOD 

>P IP  ODSINP . OLD/NV=ODS I NP . NEW ; O/RE 
>PIP  ODSINP . NEW/NV=ODSOUT . NEW » O/RE 
/RUN  ODSGRF 
>RUN  ODSDET 

>PIP  ODSIR .NEU^QDSOUT .NEW50/RE 
;RUN  CNVL 

>PIP  ODSINP . OLD/NV=QDSINP . NEW 5 O/RE 
>P I P  ODS I R . OL D/N V =ODS I R . NE W  5  0/ RE 
>P IP  ODSINP . N£W/NV=ODSOUT . NEW } O/RE 
-RUN  ODSGRF 
/ RUN  ODSIOT 

>P IP  ODSSUN . OLD/NV^GDSSUM .NEW I O/RE 
-RUN  ODSWV 

>PIP  ODSINP . OLD/NV=ODSINP . NEW J O/RE 
>PIP  ODSINP. NEU/NV=ODSOUT. NEW! O/RE 
>PIP  ODSWV. NUM/NV=ODSWV.NUM; -1/RE 
>RUN  ODSMOD 

>PIP  ODSINP. OLD/NV=ODSINP. NEWS O/RE 
>PIP  ODSINP . NEW/NV=ODSOUT . NEW ! O/RE 
>RUN  ODSGRF 
RUN  ODSDET 

>P IP  ODS IR , NEU=ODSOUT . NEW ; O/RE 
./RUN  CNVL 

>PIF  ODSINP. 01D/NV=0DSINP.NEW;0/RE 
PIP  ODSIR .  OLD/NV=ODSIR .  NEW  i  O/RE 
/PIP  ODSINP. NEW/NV=ODSOUT. NEW  5 O/RE 
/RUN  ODSGRF 
-PUN  Ol'STOT 

/PIP  ODSSUM. OLD/NV=ODSSUM. NEW! O/RE 
/ PIP  ODSINP . NEW/NV=OI*SSUM . OLD ! O/RE 
/-RUN  ODSHVY 

. PIP  ODSIR. NEW=0DS0UT.NEU;0/RE 
/RUN  CNVL 

PIP  ODSINP. 0LD/NV=0DSINP. NEW! O/RE 
-PIP  ODSIR, OLD/NV=ODSIR. NEW ;0/RE 
.PIP  ODSINP .NEW /NV=ODSOUT .NEW f O/RE 
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RUN  OUSGkF 
PUN  ODSPSP 
/  5 

»  CONTINUE  WITH  ODS?  CY/NKN 

's  f 

;;  » 

>»  >>  END  OF  ODS  RUN 

>i  » 

>5  »  (YOU  MAY  WISH  TO  PURGE  FILES  J.SAV  AND  *.PRM  AT  THIS  POINT) 

>i  » 

>* 

>0  <EOF> 

> 
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TEST  CASE  6 
Program  Users  Manual 


A  'll 


tODSRUN 


»  WARNING! 


IT  IS  NECESSARY  TO  CLEAR  ALL  ACTIVE  ODS  FILES. 


>*  OK  TO  CLEAR  FILES  AUTOMATICALLY?  [Y/N]!N 
>#  HO  YOU  WISH  TO  EXIT?  CY/N3JN 


>RUN  ODSXQT 


[[[  OPTICAL  DETECTION  SYSTEM>  V  2.0  333 
SCIENTIFIC. TECHNOLOGY  ASSOCIATES  (1981) 


> INPUT »  EDIT*  RUN?  INPUT 
>SELECT  NEW  SEQUENCE?  CY/N3!  Y 


4+4  ODS  MODULES  444 


SOURCE 

ATMOS 

OPTICS 

MODULATOR 

DETECTOR 

PREAMP 

AMPLIFIER 

PLOT 


>ENTER  MODULE  LIST  (TYPE  VEND*  TO  EXIT) 
//SOURCE .  PLOT .  PREAMP ,  PLOT .  AMPLIFIER » PLOT » /END 

/HOW  MANY  SOURCE  WAVELENGTHS  (LE  32767 ) ?  1 


...THE  MODULES  YOU  HAVE  REQUESTED  ARE! 


SOURCE 

(P) 

PREAMP 

(P) 

AMF'LIF  IER 

(P) 

>0K?  ENTER  Y  OR  N!  Y 


A-  22 


SAVE  THIS  SEQUENCE"? 

:>  u arming:  this  will  supersede  any  existing  sequence  « 

>CY/NK  Y 
>RUN  WVINP 

this  is  the  wave  form  section,  do  you  want  help?  < r/N) : 

N 

ENTER  WAVELENGTH  (IN  MICRONS)  FOR  WHICH  SIGNAL  IS  TO  BE  DETECTED: 

4. 

ENTER  NAME  OF  SIGNAL  YOU  WISH  TO  USE! 

SQUARE 

BEFORE  SELECTING  THE  TIME  STEP  TO  BE  USED.  PLEASE 
CONSIDER  THE  FOLLOWING: 

AMPLIFIER  !  ALLOWED  TIME  STEP  (DT) 


w 

1 

DT 

<= 

1.25E-04 

SEC 

H 

! 

DT 

5.00E-03 

SEC 

F 

j 

DT 

2.50E-04 

SEC 

C 

1 

DT 

<  = 

5.00E-04 

SEC 

D 

i 

DT 

<= 

1.00E-03 

SEC 

U 

! 

DT 

<= 

1.00E-02 

SEC 

NOW  PLEASE  ENTER  THE  TIME  STEP  (SEC.) 5 
0.0005 

NOW  PLEASE  ENTER  THE  UPPER  LIMIT  ON  TIME  DURATION  OF 
THE  SIGNAL  (SEC.): 

0.5 

THIS  IS  THE  SQUARE  WAVE.  PLEASE  ENTER  THE  FOLLOWING  VALUES*. 
THE  AMPLITUDE 
1.0 

THE  WIDTH  OF  THE  BASE  OF  THE  SQUARE  WAVE 

0.1 

THE  SEPARATION  BETWEEN  SQUARE  WAVE  PULSES 

0.1 

THE  LENGTH  OF  TIME  BEFORE  THE  SQUARE  WAVE  BEGINS(DEADZONE) 
0.05 

THE  BASELINE  (AMPLITUDE  MINIMUM) 

0. 

DO  YOU  WANT  TO  ADD  NOISE  TO  THIS  WAVE  (Y/N)? 

Y 

DO  YOU  WANT  [WHITE]  NOISE  OR  [1/F]  NOISE  ? 

WHITE 

YOU  WANT  NOISE  TO  BE  ADDED  TO  YOUR  CHOSEN  INPUT  SIGNAL. 

-  LEASE  FJ-ER  AN  RMS  A  .PL  I  TUBE. 

■  .A 
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-RUN  HVF'INP 

THIS  IS  THE  FRE-AMF-  SECTION.  USING  THE 
HEAVISIDE  OPERATORS.  DO  YOU  UANT  HELP  (Y/N)? 

N 

DO  YOU  WISH  TO  USE  ONE  OF  THE  SUPPLIED  ROUTINES 
(Y/N)? 

N 

DO  YOU  WANT  TO  INPUT  POLES  AND  ZEROES  DIRECTLY  (Y/N)? 

Y 

IS  THERE  A  QUADRATIC  TERM  PRESENT  (Y/N)? 

N 

PLEASE  ENTER  THE  NUMBER  OF  POLES*. 

1 

PLEASE  ENTER  THE  NUMBER  OF  ZEROES*. 

0 

PLEASE  ENTER  POLE  #1J 
490. 2 

PLEASE  ENTER  THE  VALUE  OF  A  PREFACTOR  *. 

1. 

PLEASE  ENTER  THE  UPPER  LIMIT  ON  THE  TIME  FOR  THE 
IMPULSE  RESPONSE  (E10.0).  THIS  WILL  BE  IN  SECONDS. 

SELECT  YOUR  TIME  LIMIT  TO  BE  BETWEEN  THE  FOLLOWING  TWO 
TIMES:  2.560E-01  SECONDS  AND  O.OOOE-Ol  SECONDS. 

0.02 

>RUN  HVYINP 

THIS  IS  THE  AMPLIFIER  SECTION,  USING  THE 
HEAVISIDE  OPERATORS.  DO  YOU  WANT  HELP  (Y/N)? 

N 

DO  YOU  WISH  TO  USE  ONE  OF  THE  SUPPLIED  ROUTINES 
(Y/N)? 

Y 

WHICH  RESPONSE  FUNCTION  DO  YOU  WANT?  THESE  ARE  THE  AVAILABLE  CHOICCS 
W 
H 
F 
C 
D 
G 

LP 

BP 

HP 

BR 

NOW  ENTER  YOUR  CHOICE*. 

LP 

IT  WILL  NOT  BE  NECESSARY  TO  ENTER  A  MAXIMUM  TIME. 

IT  WILL  BE  COMPUTED  BASED  ON  A  FIXED  NUMBER  OF  POINTS  =  256. 


NOW  YOU  MAY  ENTER  THE  FIRST  FREQUENCY. 

PLEASE  REMEMBER  THAT  IF  YOU  HAVE  CHOSEN  LF'  OR  HP  AMPLIFIERS! 
THAT  THIS  IS  THE  CUTOFF  FREQUENCY.  IF  YOU  HAVE  CHOSEN  THE 
BP  OR  BR  AMPLIFIERS!  THEN  THIS  IS  THE  FREQUENCY  OF  THE 
LEADING  EDGE  OF  THE  BAND.  PLEASE  ENTER  THE  1ST  FREQUENCY* 

100. 

NOW  YOU  MAY  ENTER  THE  2ND  FREQUENCY. 

PLEASE  REMEMBER  THAT  IF  YOU  HAVE  CHOSEN  LF‘  OR  HP  AMPLIFIERS! 
THAT  THIS  FREQUENCY  WITH  THE  1ST  FREQUENCY  MARK  THE 
TRANSITION  BAND.  FOR  THE  LP!  THIS  WILL  BE  >  1ST  FREQ. 

FOR  THE  HP!  THIS  WILL  BE  <  1ST  FREQ.  IF  YOU  HAVE  CHOSEN 
THE  BP  OR  BR  AMPLIFIERS!  THEN  THIS  IS  THE  FREQUENCY  OF 
THE  TRAILING  EDGE  OF  THE  BAND.  PLEASE  ENTER  THE  2ND  FREQUENCY 
120. 

PLEASE  ENTER  THE  ATTENUATION  (IN  DB)t 
30. 

!>  ! 

>*  CONTINUE  WITH  ODS?  CY/N3JY 

>  ! 

>RUN  ODSXQT 
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> INPUT !  EDIT!  RUN?  RUN 
-SELECT  NEW  SEQUENCE?  CY/N35  N 

>pip  odssum.old/nv=nl: 

>pip  odssum.neu/nv=nl: 

>PIP  ODSINF' ,OLD/NV=NL! 

>PIP  ODSINF’.  NEU/NV=NLJ 
>RUN  ODSUV 

>PIP  ODSINF’.  OLD/NV=ODSINP.  NEW  JO/RE 
>PIP  ODSINF’.  NEU/NV=ODSOUT.  NEW;  O/RE 
>PIP  ODSUV . NUM/NV=ODSWV . NUM i -1 /RE 
;  RUN  ODSGRF 
; RUN  ODSTOT 

>P IP  ODSSUM . OLD/NV=ODSSUM .  NEW } O/RE 
>PIP  ODSINF' . NEU/NV=ODSSUM . OLD! O/RE 
>  RUN  ODSHVP 

;  F  IP  ODSIR,NEU=ODSOUT. NEW? O/RE 
-  F'UN  CNVL 

/PIP  ODSINF’ . OLD/NV=ODSINP . NEW i O/RE 
/■PIP  ODSIR.OLD/NV=ODSIR. NEW! O/RE 
;PIF’  ODSI NP • NEU/NV=ODSOUT . NEW i O/RE 
-RUN  ODSGRF 


.✓RUN  ODSHVY 

>PIP  ODSIR.N£W=ODSQUT.NEWiO/RE 
>RUN  CNVL 

>P IP  ODSINP . OLD/NV=GDS I NP . NEW J O/RE 
>PIP  ODSIR.OLD/NV=ODSIR. NEW i O/RE 
>PIP  ODSINP . NEU/NV=ODSQUT . NEW » O/RE 
>RUN  ODSGRF 
>RUN  ODSDSP 

>r 

>*  CONTINUE  WITH  ODS?  CY/N35N 

>; 

>;  » 

>»  »  END  OF  ODS  RUN 

>;  » 

>;  »  (YOU  MAY  WISH  TO  PURGE  FILES  *.SAV  AND  *.PRM  AT  THIS  POINT) 

>»  » 

>; 

>8  <EOF> 

>PIP  *,SAV»*/DE 

>pip  *.prm;*/de 

> 


A  -  26 


